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FINAL REPORT 

CANDIDATE CONFIGURATION TRADE STUDY, 
STELLAR-INERTIAL MEASUREMENT SYSTEM (SIMS) 
FOR AN EARTH OBSERVATION SATELLITE (EOS) 

ABSTRACT 


A two year trade study was performed by the Charles Stark 
Draper Laboratory Division of the Massachusetts Institute of Technology 
for the NASA L. B. Johnson Space Center, under the technical direction 
of the NASA Goddard Space Flight Center. 

Ten candidate SIMS configurations were originally defined in 
the first interim technical report in November 1971 and reduced to three 
in the second such report in February 1972. The latter three were then 
studied in depth in the third interim technical report in June 1972. After- 
wards, two additional studies were performed where one involved a more 
detailed error analysis of the SIMS-A candidate and the other was an 
investigation of the use of known landmarks in the payload sensor imagery 
to estimate spacecraft attitude in the event of star sensor failure for the 
SIMS-A and SIMS-B candidates. The results of these two studies were 
reported on in the Final Report in January 1973. 

This addendum to the Final Report gives the results of the latest 
study which was an investigation of the effects of spacecraft orbital 
ephemeris and attitude errors on the ability to determine the locations of 
unknown landmarks with respect to known landmarks in the payload sensor 
imagery. 


by Robert White 
Fred Grant 
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PREFACE 


This addendum to the Final Report presents an investigation 
of the effects of spacecraft orbital ephemeris and attitude errors on the 
ability to determine the locations of unknown landmarks with respect to 
known landmarks in the payload sensor imagery of an Earth Observation 
Satellite in a sun -synchronous orbit at 540 nautical miles altitude. Most 
of the effort was concerned with the ability to determine landmark lo- 
cations in an observational pass over the continental USA. 

The data in this report gives the sensitivity of landmark lo- 
cation error to errors in spacecraft ephemeris, attitude, gyro bias drift, 
and higher harmonics in the Earth's gravitational field. Also presented 
is data showing the relative roles that the several error sources have on 
landmark location determination. Finally, some results are given show- 
ing the effects of typical sets of errors in the smoothed estimates of 
attitude and gyro bias drift for SIMS-A. It has been found that the errors 
in the smoothed estimates of attitude and gyro bias drift do not affect 
landmark location determination as much as some people might suspect. 
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SECTION 1 


INTRODUCTION AND SUMMARY 
1.1 INTRODUCTION 

This report has been prepared as an addendum to the previously 
published Final Report and covers the studies performed between 
1 April 1973 and 30 June 1973 by the Charles Stark Draper Laboratory 
Division of the Massachusetts Institute of Technology (MIT/CSDL) on the 
"Candidate Configuration Trade Study — Stellar-Inertial Measurement 
System (SIMS) for a Proposed Earth Observation Satellite (EOS) 1 ' for the 
NASA Goddard Space Flight Center (GSFC). 

Three prior Interim Technical Reports 85, 141,169 and a Fi n a.i 
177 

Report have been published. Excerpts from MIT/CSDL Technical 
Proposal No. .71-173, dated June 1971, including the statement of work 
for the first eleven months of this effort, were provided as Appendix A of 
Reference 85. , Excerpts from MIT/CSDL Technical Proposal No. 72-176, 
dated 16 May 1972, including the statement of work for the following ten 
months of the effort, were provided as Appendix E of Reference 169. The 
final three months of this effort were spent on additional studies proposed 
in MIT/CSDL Technical Proposal No. 73-132 and reported on in this 
addendum to the Final Report. Two new studies were to be performed in 
this final three month effort which are as follows: 

1) Landmark Location Study - Investigate the effects of space- 
craft orbital error, in addition to attitude determination 
errors, on the ability to determine the locations of unknown 
landmarks with respect to known landmarks in the payload 
sensor imagery. 


* Superscripts refer to similarly numbered references in the reference 
section (Section 6) of this report or the reference section of the previous 
reports of the SIMS Trade Study. 



2) Active and Passive Landmark Mechanization Study - Perform 
a preliminary survey of the various techniques or devices on 
ground which might be suitable for obtaining measurement 
data for EOS orbit and attitude determination. 

The data presented in this addendum to the Final Report re- 
presents the results of the Landmark Location Study. Due to other 
commitments of key personnel, MIT/CSDL was unable to perform the 
Active and Passive Landmark Mechanization Study during this three month 
period. Separate negotiations are now underway to perform this study 
at a later time. 


1.2 BACKGROUND OF SIMS TRADE STUDY 

Four categories of SIMS candidate configurations were originally 
required to be evaluated and compared in the SIMS Trade Study. These 
were designated as A, B, C, and D. Later, E was added and D was sub- 
divided into sub-groups. Thus at the time of the First Interim Technical 
Report, 10 candidates were defined as potential SIMS design approaches. 
These were eventually narrowed down to 3 candidates in the Third Interim 
Technical Report and were simply designated as SIMS-A, SIMS-B, and 
SIMS-D. The accuracies given for the three candidates in the Third Interim 
Technical Report were obtained by performing a covariance analysis with 
an optimal smoother. This data provided a statistical indication of the 
accuracy possible after smoothing star and gyro rate measurements over 
a number of satellite orbits. In the error studies presented in the Final 
Report a more detailed error study was performed on the SIMS-A can- 
didate. Smoothed results were obtained for the estimates of spacecraft 
attitude and gyro bias drift and also for the uncertainties in these estimates. 
The Final Report also presented the results of an investigation into the 
use of known landmarks in the payload sensory imagery to estimate 
attitude in the event of star sensor failure in SIMS-A and SIMS-B. 
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1.3 SUMMARY 

This subsection briefly summarizes the principal data 
given in this report. This report was primarily concerned with 
determining how errors in the estimation of spacecraft attitude and 
orbital ephemeris would affect the ability to determine the locations of 
landmarks in the payload sensor imagery during an observational pass 
over the continental USA. The effect of incomplete modeling of the 
Earth’s non -spherical gravitational field and of the uncertainties in the 
harmonic coefficients of that field are also studied. The satellite orbit 
was assumed to be a sun- synchronous orbit with an altitude of 540 nautical 
miles and an inclination of 99 degrees. The scan beam of the payload 
sensor generated a swath of imagery 100 nautical miles wide by sweeping 
back and forth across the satellite ground track to a maximum angle of 
5.3° to each side of the orbital plane. 

Section 2 of this report presents the analytical and supporting 
equations of this study. Section 3 gives the error study results which 
are roughly divided as follows: 

1) Sensitivity Studies - Data is given showing the sensitivity 
of landmark error to error in each component of satellite 
attitude, orbital ephemeris, and gyro bias drift. Included 

also are the effects of attitude libration and various gyro errors 
and noises. 

2) Gravitational Modeling - Data is given showing the landmark 
errors caused by the higher harmonics in the Earth's gravi- 
tational field. Included also are the effects of uncertainties 
in these harmonics. 

3) Effect of Smoothed Estimates of Attitude and Gyro Bias Drift - 
Data is given showing the landmark errors caused by errors 
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in the smoothed estimates of attitude and gyro bias drift 
of previous SIMS studies. 

4) Effect of Combined Errors in Spacecraft Attitude and Orbital 
Ephemeris - Data is given showing the combined effect of 
nominal errors in satellite attitude and orbital ephemeris. 

Section 4 presents some of the conclusions drawn from this study 

and Section 5 gives some recommended future studies. In Appendix A 

177 

some corrections are given for the previously published Final Report 
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SECTION 2 


LANDMARK LOCATION STUDY ANALYSIS 
2.1 INTRODUCTION 

The primary objective of the landmark location study was to de- 
termine the effects of spacecraft orbital ephemeris and attitude errors on 
the ability to determine the locations of unknown landmarks with respect to 
known landmarks in the payload sensor imagery. Sufficient data was gen- 
erated to show the sensitivity of landmark location determination to errors 
in each component of the spacecraft orbital ephemeris and attitude. This 
data""can'be of greafvalue in indicating what accuracies“are r equi red”irTtfie 
orbital ephemeris and spacecraft attitude in order to achieve a given 
accuracy in landmark location determination. No attempt was made to 
correct (or update) our knowledge of the orbital ephemeris and spacecraft 
attitude using the observations made on known landmarks since this is a 
different problem and should be given separate treatment in future studies 
using optimal filter and/or smoother techniques (see Section 5). 

In the present study the observed discrepancy in the location of 
a known landmark is used to correct our knowledge of the locations of 
other arbitrary but unknown features seen at different times in the pay- 
load sensor imagery. This gives us an indication of how well we can de- 
termine the location of an unknown landmark (or feature) with respect to 
a known landmark. During the period between observations of the known 
and unknown landmarks, the error in each component of the satellite orbital 
ephemeris and attitude will change in a prescribed manner; and it should be 
noted that it is the 'change' in the errors of spacecraft orbital position and 
attitude, and not the errors in these quantities, which will have the most 
significant effect on the relative landmark location accuracy. 

One other objective of the landmark location study was to determine 
the effects of uncertainties in the Earth's gravitational model on the 
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propagation of the orbital ephemeris for an Earth Observation Satellite 
(EOS). Coupled with this was a desire to determine how complete the 
gravitational model should be in the landmark location problem. Today, the 
Earth's gravitational model is fairly well known and can be represented by 
a large number of terms in a series expansion. Most of these terms are 
relatively small and are not required in this particular application. 

In this study most of the emphasis was placed on landmark location 
accuracy within the continental USA (between 30° and 50° latitude). Figure 
2-1 shows a typical observational pass over the USA, The EOS satellite is 
assumed to be in a circular orbit with an altitude of 540 nautical miles and 
an inclination of 99 degrees. Onboard the satellite there is assumed to be a 
thematic mapper (or multi -spectral scanner) whose beam is directed down- 
wards and sweeps back and forth across the ground track to generate a 
swath of imagery 100 nautical miles wide. Under nominal attitude conditions, 
the spacecraft body axes (Xg, Yg, Zg) are as shown in Figure 2-1, with 
Xg in the direction of satellite motion and Z R pointing to nadir. The 

scanner beam sweeps back and forth about X D with a maximum excursion 
o 

of 5. 3 with respect to Zg. The angle between the beam and Z R is denoted 
as e and is positive when in the same direction as Y R . 

Most of the performance results of this study are for an observa- 
tional pass which starts at 50° latitude with specified errors in satellite 
position, velocity, attitude, and gyro bias drift. The satellite position and 
velocity at any later point are obtained by solving the equations of motion. 

The spacecraft attitude at any later point is obtained by using a fourth order 
Runge Kutta integration of the indicated body rates of three body mounted 
gyros. The input axes of the gyros are assumed to be oriented in the same 
way as in SIMS-A and SIMS-B, which is parallel to the spacecraft body 
axes Xg, Y r and Z R . A realistic simulation was used for the gyros, in- 
cluding the effects of bias drift error, random drift, scale factor error, 
input axis misalignment, and quantization. 

It is important to note that no updates were made in the satellite 
orbital ephemeris and attitude during the relatively short observational pass 
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over the USA since it was felt that such updates are usually random in 
character and would probably introduce more error in the relative land- 
mark location determination than the changes in the errors of the assumed 
orbital ephemeris and attitude. A possible exception to the above might be 
the case where use is made of a smoothed well-behaved satellite position 
and attitude history based on batch processing. 

It should also be noted that no consideration was given to the effect 
of thematic mapper errors on landmark location determination since this was 
beyond the scope of this effort. 

2. 2 RELATIONSHIP BETWEEN LANDMARK LOCATION AND 

SATELLITE ATTITUDE AND POSITION 

It is desired to obtain the equations relating spacecraft attitude 
and position to landmark location. We first obtain equations for the land- 
mark position vector as a function of the orbital position vector and the line 
of sight vector of the scan beam at the time of observation of the landmark. 

Figure 2-2 shows the beam vector, s, the landmark position vector, 
J., and the satellite position vector, r. Also shown are the satellite altitude, 
H, and the instantaneous deviation angle, 5 , of the scan beam vector from 
the satellite position vector (or local vertical). From the figure we have: 

]_=£+£ ( 2 - 1 ) 

where all vectors are assumed to be in the basic inertial coordinate system, 
which is Earth centered with its Z-axis along the north pole and the other 
two axes are in the equatorial plane. 

Since the deviation angle, 8, will be small in the present application, 
a flat Earth can be assumed and Equation 2-1 can be expressed as follows: 

2 = r + (H/cos8)u LOS (2-2) 
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SATELLITE 




where 


—LOS = UNIT(s) 

cos 6 = u lJDS • [ - UNIT ( r ) ] 

The direction of the scan beam in spacecraft body coordinates 
at any given instant is completely defined by the scan beam angle, t , as 
shown in Figure 2-3. In basic inertial coordinates, the unit vector defining 
the direction of the beam is given by: 

u los = Eb cos e + Zb sin e (2-3) 

where ^ and are unit vectors defining the directions of the Y and Z 
axes of the spacecraft body in basic inertial coordinates and are obtained 
as follows: 



0 


o~ 

£b t io t ob 

1 

0 

* -B ~ T IO T OB 

0 

1 


where T OB is the transformation matrix from body to orbital coordinates 
and Tjo is the transformation matrix from orbital to basic inertial co- 
ordinates. The relationships between the Body, Orbital and Basic Inertial 
Coordinate Systems are shown in Figures 2-4 and 2-5. The equations for 

T IO and t ob are: 


cn 

-SO ci 

sn Si 

sn 

cn Ci 

-cn Si 

o 

si 

ci 


(2-5) 
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(North Pole) 



Figure 2-4 Basic Inertial and Orbital Coordinate Systems 



Figure 2-5 Orbital and Spacecraft Body-Fixed 
Coordinate Systems 
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-C0 CO 


-ci jj s o ~stfj so 

+ s\p S<f) CQ -Clp S<p CO 


OB 


cijj co sip co ~c4> so 

+ s p s <t> so -cijj s<j> so 


Sp C(f) 


-op cp 


sp 


( 2 - 6 ) 


where s and c denote sine and cosine; flis orbital right ascension; i is 
orbital inclination ; and Q, p, and i p are the spacecraft pitch, roll and yaw 
angles, respectively. Substitution of Equation 2-6 into Equation 2-4 
gives: 


Xb 



-sip so - cip sp co 
sip c$ - a p s p so 
- cip cp 




-C(f) co 

-c p so 

S 0 


(2-7) 


( 2 - 8 ) 


Note that Equations 2-7 and 2-8 may be used to compute either the 
true or estimated body axis vectors depending on whether the true or 
estimated spacecraft attitude angles are used. 


When the expression for u^g in Equation 2-3 is substituted into 
Equation 2-2, the following is obtained for the landmark position vector: 


_1 = r_ + (H/cos6) Zg cose + ^ sin ej (2-9) 

where j. can either be the true or estimated landmark position vector depend- 
ing on whether true or estimated values are used for r, z^, and y,-,. In the 

— — -D — JJ 

present effort, both values of_l are computed for each landmark. The dif- 
ference between the two values represents the error in estimating the 
landmark location due to errors in the estimates of orbital position, r, and 
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spacecraft attitude, z^ and y^. In this study, most of the attention has 
been directed towards the horizontal components of the landmark error. 
These components may be computed relative to north-south and east-west 
axes from which the errors in latitude and longitude may be computed. How- 
ever, the performance results in this report have been presented in cross - 
track (normal to the orbital plane) and down-range coordinates. In the 
tabular results these two terms have been shortened to 'track' and 'range', 
respectively. The track and range components of the landmark error were 
chosen in order to clarify certain interesting effects which are discussed 
later in the section on results (Section 3). The track and range components 
of the estimated or true landmark position vector may be computed in true 

[UNIT (vx r) 

[UNIT (v) ] 

where r and v are the true satellite position and velocity vectors, respect- 
ively, and l_ is the estimated or true landmark position vector. In this case 
the track and range components of the error in the estimate of landmark 
location are: 

e TRK = 1 TRK E 

e RNG = Wg e 

where the subscripts E and T denote the estimated and true values, 
respectively. 

In Section 2. 6 a discussion is given of the technique used to 
correct the position estimates of arbitrary features (or unknown landmarks) 
observed during the observational pass. This correction is based on the 
values of and e^^, obtained for a single control point observed 

during the pass. After correction, the errors in the position estimates 
of the arbitrary features are considered to be an indication of how well 
these features can be mapped with respect to the control point. 


TRK r 


TRK r 


( 2 - 12 ) 

(2-13) 


( 2 - 10 ) 

( 2 - 11 ) 


local vertical coordinates as follows: 


X TRK “ R E sin 1 


X RNG = R E sm 


UNIT (1) | * 
UNIT (1 ) 
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The latitude and longitude of the landmark vector, 1, can be com- 
puted as follows: 


LAT 


tan 


-1 


Y + i 2 


LONG = SIGN (1 Y ) cos 


-1 


X 


l 2 + l 2 
X ' l Y 


(2-14) 

(2-15) 


where the subscripts denote the components of 1 in Earth fixed coordinates, 


2. 3 PROPAGATION OF STATE ESTIMATES OF SPACECRAFT 
ATTITUDE AND GYRO BIAS DRIFT 

1 89 

In the Third Interim Technical Report and in the Final 
177 

Report , considerable data was presented on the uncertainties and 
errors in the smoothed estimates of spacecraft pitch/ roll, and yaw, and 
bias drift rate of each gyro. These errors were for a particular specified 
time in the data processing interval, usually the mid-point of the data 
interval. 


The ground-rule adopted in the present error study was that there 
be no update of spacecraft attitude and gyro bias drift during the observational 
pass over the USA. At the start of an observational pass a given set of 
values (including errors) was assumed for pitch, roll, yaw, and gyro bias 
drift. In some cases these values represented the smoothed results of pre- 
vious SIMS studies. Thereafter, during the observational pass, the space- 
craft attitude (6, <ft and i p) was obtained by performing a fourth order 
Runge Kutta integration of the following rates: 


e 


Sl///c</> 

-Cl/// C(f> 

0 

</> 

= 

Cl// 

si ft 

0 

A 


S<f> Si// 1 C(f) 

-S</> Cl/z/ C(f) 

-1 


-y 


(2-16) 
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where o> M is the vector for the measured gyro rates which includes the 

effects of true gyro bias drift, bias and random errors in gyro scale factor, 

errors in gyro input axis alignment, and errors due to random drift (see 

177 

subsection 2. 2. 2 of the Final Report ). The vector b represents the 
estimate of bias drift for each gyro, which is used to compensate for the 
bias drift present in 

It is important to note that the spacecraft attitude was propagated 
in the same manner as was done in the detailed error analysis of SIMS-A 
except that no updates were made with star measurements. 

The true spacecraft attitude history used to generate most of the 
results in the present effort was nominal (i.e., the roll and yaw angles 
(rf) and *]/) were zero, and the pitch angle (0) changed at a constant orbital 
rate). However, in some cases, use was made of the following attitude 
history: 

= <*>t + 9^ sin ( t + 120°) + 

(f > T = <t> L sin ( cu L t) + 4> 0 (2-17) 

^T = sin * W L * + 240 °) + 

where 0^, and ^ are the true pitch, roll, and yaw angles, re- 

spectively; a; is the constant orbital rate; t is the time since t = 0; 0 q, 

</> q , and if/q are the attitude angles at t = 0 ; 0 ^, and ^ are the 

libration amplitudes (0.6° each); and is the libration angular rate (0,478 
degrees per second). 

2. 4 PROPAGATION OF ORBITAL EPHEMERIS 

An important, if not most significant, source of error in the 
present landmark location problem is the error in the knowledge of 
satellite position and velocity. This error is primarily due to errors in 
tracking system instrumentation and uncertainties in gravitational model- 
ing, atmospheric drag, and tracking station location. 
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Various values have been quoted for the uncertainties in orbital 
ephemeris depending upon the extent and quality of ground tracking and 
data processing. For the purposes of this study it was indicated by GSFC 
that a satellite position uncertainty of 50 meters should be considered as 
nominal with most of this being in the down-range component. Based upon 
previous CSDL experience with uncertainties in orbital position and 
velocity, the following uncertainties were selected as being nominal for 
the components of position and velocity; 50 meters in down-range; 20 
meters in altitude and cross-track; 0. 05 meters per second in altitude 
rate; and 0.02 meters per second in cross-track and down-range rate. 
Actually, the selection of the above values as nominal is not too important 
since sufficient sensitivity results were generated to show the effect of 
different values of error in each component of the orbital ephemeris. 

In the present study, a given set of errors in orbital position and 
velocity were assumed at the start of an observational pass over the USA. 
These errors were then propagated along the orbit to the points of interest. 

There are several methods by which the errors in position and 
velocity may be computed for a later point in orbit. The one employed in 
the present study used a fourth -order Runge Kutta integration of the follow- 
ing equations of motion to obtain the position, r, and velocity, v, at some 
desired point for the case of no initial errors and for the case with initial 
errors. The difference between the two solutions represents the error at 
that point. 

• r_ 

v = - /x (2-18) 

r = v (2-19) 

where ^ is the gravitational constant. 

Actually, there are simpler ways to compute the position and 
velocity at some later point in orbit when only a central force field is 
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assumed for the Earth. One of these is the following closed-form vector 
solution: 


r = 


v = 


1 - ^ [1 - COS (A0)| 


r„ + r r 0 sin (A0 ) 

“° JTT* 


p r i 


0 [ 1 - cos (A0)] - p y p sin (A0 ) 

(A0)] 


—0 


+ 


1 - — 1 1 - cos 


^0 


( 2 - 20 ) 


( 2 - 21 ) 


where r^ and are the initial position and velocity vectors; r and v are 
the position and velocity vectors after a true anomaly change of A0; and 
p is the parameter of the orbit. Although this method was not employed in 
the present effort, it has been presented as a point of interest. 


A third method of computing the errors at some later point is to 
use the following approximate error propagation equations which are 
satisfactory for a circular or near circular orbit. In addition to saving 
computational time, the following equations are useful in demonstrating 
certain relationships between the errors: 


"ALT 


e ALT 0 (2 - coswt) + e RNG 0 SinWt 


+ T (1/co) sin CD t + e RNG (2/co)(l - coscot) (2-22) 


RNG 


e ALT sinojt ’ 3wt ) + e RNG (2 cos cot - 1) 


e ALT (2/co) (1 - cos cot)+ e RNG (1/co) (4 sin cot - 3 cot) 

0 ‘ 0 (2-23) 


e TRK = e TRKg cos cot + e TRK ^ (!/co) sin cot 


(2-24) 
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where the subscripts ALT, RNG, and TRK denote the altitude, down-range, 
and cross-track components of the error, respectively; the subscript 0 
denotes the initial error at t = 0; and co is the spacecraft orbital rate. 

Note that Equations 2-22 through 2-24 give only the position errors at a 
later point in orbit since the velocity errors at that point do not affect the 
landmark location determination at that point. Later, in the discussions 
of Section 3, it will be found that these equations are very useful in ex- 
plaining the performance results. 

2. 5 CONSIDERATION OF EARTH'S GRAVITATIONAL MODEL 

In Section 2.4 consideration was given to propagation of the 
orbital ephemeris errors in the Earth's central force field (i. e. only the 
gravitational constant v was assumed to be present). Although most of 
the performance results of this study were for this case, some data was 
generated showing the effects of higher harmonics in the Earth's gravita- 
tional potential. The primary reason for considering a more exact model 
of the Earth's gravitational potential was to determine how complete the 
model should be for the particular problem of landmark location deter- 
mination in the continental USA. 


The following material presents some of the basic equations used 
for modeling the Earth's gravitational field. These equations are based 
upon some important work by S. Pines in Reference 178. 


The gravitational potential, U, can be expressed as: 

oo n 


u = JL 


1 +'S 'S (u) [c„ .cos niA + S sin mxl 

Z-, \ r / n, m [ n, m n, m J 


n= 2 m=0 


(2-25) 
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where 


or 


C and S are constants 
n, m n, m 




= 0 ) 


P n,m< u) 


(1 - u 2 ) /dT +m 
2 n n ! \diy 



u = z/r = cos (Latitude = 90° - (j>) 

X = tan” 1 (y/x) = Longitude 
a = Equatorial radius 
r = Radial distance 
H = Gravitational constant 

The general equation for the accelerating force vector is 
f = dU 

— dx (2-26) 


aU dr , aUdH. 4 - c*J d>L 
— " 5r dx ^u dx <9A dx 


(2-27) 


Because of the existence of a singularity at sin <f> = 0 that develops 
with the use of these equations, it was found convenient to change to a four 
parameter system. For this we have 


r = 


(2-28) 


where s, t, and u are the direction cosines of r. The following general 
expression for the force vector can then be obtained: 


f = ~ — ^ iir + f a r -r + a l — + a 2 1- + a 3 - 1 (2-29) 

r 1 J 

- conic acceleration + disturbing acceleration 


2-15 



After evaluating the terms of the disturbing acceleration and 
reducing, the following expression for a is obtained: 


co n 


O' 


S S n+2 A n+1, m+1 R m C n, m + r m S n, m ] 


( 2-30 


n= 2 m=0 


where 


R 


m 


/, 2, m/2 , 

= (1 - u ) cos mX 


I m = (1 - u 2 ) mf2 sin mX 


( 2 - 31 ) 

( 2 - 32 ) 


and where a recursion equation for A is 

^ n, m 


A 


n+1, m+1 uA n, m+1 + ( n+Tn+1 ) A n# 


m 


( 2-33 


The equations for the other o' s (i.e. and « 3 ) are 

quite similar in form, varying principally in the subscripts used for A. 
R, and I, 


An inspection of the above equations for the disturbing acceleration 
shows that the accuracy of the earth's non-spherical gravity model depends 
upon the number of C and S coefficient pairs that are incorporated. Likewise 
it is clear that uncertainties in the values used for the individual C and S 
coefficients will be reflected as uncertainties in the earth's disturbing 
acceleration due to the earth's non-spherical and non-symmetrical charac- 
teristics. 


In the next section, data on uncertainties in these coefficients is 
presented. Also shown is the way in which these coefficients affect the 
propagation of the orbital ephemeris as reflected in landmark errors. In 
these studies it was assumed for convenience that the initial ephemeris 
errors were zero at the start of each observational pass over the USA. 
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2.6 


LANDMARK ERROR CORRECTION TECHNIQUE 


In the present study a relatively simple correction was applied 
to the position estimates of unknown landmarks observed during an 
observational pass over the USA. This correction was assumed to be that 
which is required to offset the discrepancy found between the estimated 
and true positions of a single control point (or known landmark) observed 
during the pass, usually at the start of the pass on the northern USA border. 
The discrepancy in the position estimate of the control point was due en- 
tirely to errors in the knowledge of satellite orbital position and attitude 
at the time of observation of the control point. If it were possible to keep 
the orbital position and attitude errors fixed during the observational pass, 
then the correction would be essentially the right one for all of the features 
observed during the pass. However, these errors do vary with time and 
some cause the correction to become less effective as the distance is in- 
creased between the control point and the other features. The correction 
applied to each unknown feature during an observational pass was that 
required to correct for the errors given for the control point in Equations 
2-12 and 2-13. 

It should be noted that all of the performance results in this re- 
port are based upon the assumption of only one control point per pass. By 
using additional control points it is obvious that a better landmark correction 
could be obtained. However, the present study was limited to a single control 
point for reasons of simplicity, and shows what the performance would be 
for this rather extreme case. In the performance results of Section 3, one 
can see some of the advantages to be gained by using additional well placed 
control points. However, no real consideration was given to multiple 
control points since this was considered to be more appropriate for future 
studies. 
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SECTION 3 


PERFORMANCE RESULTS OF LANDMARK 
LOCATION STUDY 


3.1 INTRODUCTION 

The performance results in this report were primarily concern- 
ed with, first, the sensitivity of landmark location error to errors in 
each individual component of the orbital ephemeris and satellite attitude, 
and second, the combined effect of selected nominal error values for the 
components of the orbital ephemeris and satellite attitude. In addition, 
performance results are given showing the effect of incomplete gravi- 
tational modeling and of uncertainties in the gravitational coefficients. 
Performance results were also generated for various sets of errors 
obtained in the smoothed estimates of attitude of previous SIMS studies. 

Most of the results were generated for an observational pass over 
the continental USA starting at 50° north latitude, 95° west longitude, and 
ending at 30° north latitude (see Figure 2-1). The time associated with 
the above latitude change was about 358 seconds. The orbital altitude was 
540 nautical miles and the inclination was 99°. One exception to the mid- 
USA passes was for the gravitational modeling studies where the effect of 
variation in the initial longitude was investigated and found to be negligible 
for the intended application. 

Unless otherwise stated, the true spacecraft attitude history was 
one where the roll and yaw angles were zero and the pitch angle changed 
at a constant orbital rate. In those cases where attitude librations (i. e. 
oscillations) were present, the attitude history was that given in Equation 
2-17 of Section 2. 3. 
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Unless specifically stated, the gyros were assumed to be 
error-free. When gyros errors were present, the nominal values used 
for each gyro were as follows: 


True Bias Drift 
Scale Factor Bias 
Input Axis Misalignment 
(About each axis normal to gyro 
input axis) 

Random Drift (white noise (la) ) 
Quantization (la) 


0.03 deg/hr 
10 PPM 


10 arcsec 
0. 01 deg/hr 
0. 1 arcsec 


When gyro noises (i. e. random drift and quantization) were present, 
these noises were generated on a Monte Carlo basis. RMS data was 
obtained for only 20 computer runs and should therefore be considered 
as being representative. 

Most of the results in this section give the data on landmark 
error relative to the ground track. The down-range component represents 
that component parallel to the ground track, and the cross-track com- 
ponent is that component normal to the ground track. (For brevity, these 
components are referred to as 'range' and 'track' in the tables. ) 

3. 2 LANDMARK ERROR SENSITIVITY RESULTS 


The performance results showing the sensitivity of landmark 
location error to error in each component of the orbital ephemeris and 
satellite attitude are given in Tables 3-1 through 3-17. Tables 3-1 through 
3-9 give the sensitivity results for initial errors in spacecraft attitude 
and gyro bias drift, and Tables 3-10 through 3-17 give the sensitivity 
results for initial errors in the orbital ephemeris. In addition to the 
above sensitivity data, these tables also include data on the effects of: 
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1) Change in latitude 

2) Change in scan beam angle (e ) 

3) Spacecraft attitude librations 

4) Gyro errors and noise 

All of the sensitivity tables give the total (i. e. uncorrected) 
error in the determination of landmark position for a landmark observed 
at the start of the observational pass and for one or more landmarks 
observed later in the pass. * With the exception of Tables 3-7 and 3-16, 
which consider the effect of variation in the scan beam angle, e, all of 
the landmarks were observed when e - 0 ( i. e. downwards along the local 
vertical except for the cases with attitude libration in Tables 3-8 and 
3-17), In all of the sensitivity results the location of the landmark 
observed at the start of the pass (i. e. at 50° latitude) was assumed to 
be known and was determined by using perfect knowledge of the orbital 
ephemeris and satellite attitude. The discrepancy between the known 
location and that computed with imperfect knowledge of orbital ephemeris 
and satellite attitude represents the landmark error given for that land- 
mark in the tables. Using this discrepancy in the technique described in 
Section 2. 6, a correction was generated for each landmark observed later 
in the observational pass. In the tables both the uncorrected and corrected 
errors are given for these later landmarks. The corrected data gives an 
indication of how well one can map one point with respect to another in 
the imagery. 

In this report the known landmark is often referred to as the 
’control point' in order to distinguish it from the other landmarks whose 
locations we wish to determine. 


* It should be noted that all landmarks in this study were generated 
artificially as a matter of convenience. 



3. 2. 1 


LANDMARK ERROR SENSITIVITY TO INITIAL ERRORS IN 
SPACECRAFT ATTITUDE AND GYRO BIAS DRIFT 


Tables 3-1 through 3-6 give the sensitivity of landmark errors 
to initial errors in pitch, roll, yaw, and the bias drift of each gyro. All 
of these tables give data for three error magnitudes in order to show 
how the landmark error varies with each error source. A fairly good 
check can be made on the data in Tables 3-1 through 3-6 by using the 
following equations for a near nominal attitude history: 


* 0 2 = 
8 L = 


80. + At 5b 

1 y 

c(A s(A0)8 « jj 
s(A0)5^ 1 + C (A0)8^ 1 


s(A0) 

5b - 

1 - c(A0)' 

CO 

X 

CO 

[l - c 


, fs(A0) 


CO 


X 


CO 


(3-1) 

5b 


(3-2) 

5b z 

(3-3) 


where 80^ , 8(f> fib, 5 b , and 5 b are the initial errors (at 50° 

latitude) in pitch, roll, yaw, and the bias drift of the roll, pitch, and 
yaw gyros, respectively; 50 2 ,5</> 2 , and 5^ 2 are the pitch, roll, and yaw 
errors at some later point in orbit; At is the time lapse between the two 
points in orbit; co is the orbital rate; and A0 is the angular change in 
orbital position (i. e. true anomaly). At a given point in orbit, the pitch, 
roll, and yaw errors (50 2 ,5<£ 2 , 5^) will cause the following landmark 
errors for a landmark observed near local vertical during a nominal 
attitude history: 


Down-Range Error 
Cross-Track Error 



(3-4) 

(3-5) 
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where K is the distance subtended per unit angle (15.9 feet per arcsec 
at 540 nautical miles). 

In the following subsections, frequent reference will be made to 
the above equations when discussing the effects of each error source. It 
should be noted that a quick comparison between the effects of the different 
error sources can be obtained by using Table 3-7. 

3 . 2 . 1 . 1 Sensitivity to Pitch Error 

In Table 3-1 it is seen that the landmark error caused by a given 
initial pitch error is independent of latitude change. This is due to the 
fact that the pitch error remains fixed throughout the orbit and causes a 
constant angular deviation of the downward pointing Z B axis (see Figure 

2- 1) from the vertical in the orbital plane. Hence the resulting landmark 
error in range must remain constant. On the right side of Table 3-1 it is 
seen that the landmark error after correction for the initial error (using 
the known landmark at 50° latitude) is zero. This indicates that we can 
do a very good job of mapping one point with respect to another in the 
presence of pitch error. 

3.2.1. 2 Sensitivity to Roll Error 

In Table 3-2 the effect of roll error on landmark error is ini- 
tially the same in magnitude as that due to pitch error except that it is 
in track (see Equation 3-5). However, as the satellite moves ahead, the 
roll error changes in accordance with the first term on the right of Equation 

3- 2. The landmark error after correction for initial error increases 
approximately as (1 - cos A0). In Table 3-2 it is seen that for small 
latitude changes, such as across the USA, the landmark error after 
correction is very small in comparison to the uncorrected error. 
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TABLE 3-1 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL PITCH ATTITUDE ERROR 


Initial 

Final 

Landmark Error (feet) 

Landmark 

Pitch 

Latitude 

At Initial 

At Indicated 

Error After 

Error 



Latitude (50°) 

Final Latitude 

Correction (feet) 

(arcsec) 

(deg) 

Track 

Range 

Track 

Range 

Track 

Range 

1 

45 

0 


-15. 9 

■ 


-15. 9 

0 

0 

10 





-159 



-159 





50 





-795 



-795 





1 

40 



-15. 9 



-15. 9 





10 





-159 



-159 





50 

— 





-795 



-795 





1 

30 



-15. 9 








10 





-159 








50 

( 




-795 








1 

20 



-15. 9 



-15. 9 





10 





-159 



-159 





50 





-795 



-795 





1 

1 

0 



-15. 9 



-15. 9 





10 





-159 



-159 





50 





-795 



-795 





1 

0 




-15. 9 



-15. 9 





10 





-159 



-159 





50 

1 



t 

-795 



-7 95 






Notes: 1) e = 0 

2) No spacecraft librations 
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TABLE 3-2 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL ROLL ATTITUDE ERROR 


Initial 

Roll 

Error 

(arcsec) 

Final 

Latitude 

(deg) 

Landmark Error (feet) 

Landmark 
Error After 
Correction (feet) 

At Initial 
Latitude (50°) 

At Indicated 
Final Latitude 

Track 

Range 

Track 

Range 

Track 

Range 

1 

4 

5 

-15. 9 


0 

-15. 8 

C 


0 . 1 



10 



-159 



-158 



1 



50 



-795 



-792 



3 



1 

40 

-15. 9 






0. 3 



10 



-159 






2 



50 



-795 



mm 



13 



1 

30 

-15.9 



-14. 9 



1.0 



10 



-159 



-149 



10 



50 



-795 



-745 



50 



1 

2 

0 

-15. 9 



-13. 7 



2.2 



10 



-159 



-137 



22 



50 


1 

-795 



-685 



110 



1 

10 

-15. 9 



-12. 0 



3. 9 



10 



-159 



-120 



39 



50 

1 


-795 



-600 



195 



1 

0 


-15.9 



-10.0 



5. 9 



10 



-159 



-100 



59 



50 

1 


-795 



-500 

. 1 


295 

1 

r 


Notes- 1) € = 0 

2) No spacecraft librations 
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3. 2. 1.3 Sensitivity to Yaw Error 


In Table 3-3 it is seen that an initial yaw error does not pro- 
duce any error in the initial landmark (or control point) since the land- 
mark is directly below the satellite (i. e. e =■ 0). However, as the satellite 
moves ahead, the initial error in yaw causes an error in roll as indicated 
by the second term on the right of Equation 3-2. The reason for this 
behavior is that the initial yaw error should really be looked upon as an 
error in inertial space. As the satellite's roll and yaw axes rotate in 
inertial space, this error is resolved into roll and yaw components as 
shown in Equations 3-2 and 3-3. A yaw error at one point in orbit becomes 
entirely a roll error after 90° of orbital motion. In Table 3-3 it is seen 
that for landmarks observed later in the pass, the error after correction 
is the same as before correction since no correction w r as called for in 
the control point. 

In comparing the results of Tables 3-1 through 3-3, it is seen 
that yaw error has a much greater effect than pitch or roll error on the 
relative mapping capability. 

3. 2. 1.4 Sensitivity to Gyro Bias Drift Error 

The sensitivity of landmark errors to initial errors in the know- 
ledge ot bias drift of the roll (x), pitch (y), and yaw (y) gyros are shown 
in Tables 3-4 through 3-6. In each case there was no landmark error at 
the start of the observational pass since there was no error in spacecraft 
attitude at that point. However, as the satellite proceeds ahead, attitude 
errors will be introduced by the initial bias drift errors as shown in 
Equations 3-1 through 3-3. Since there was no error in the control point, 
there is no error correction for landmarks observed later in the pass. 


In Table 3-4 the results are given for an initial error in roll (x) 
gyi o bias drift. Here it is seen that the landmark error in track increases 
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TABLE 3-3 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL YAW ATTITUDE ERROR 


Initial 

Final 

Landmark Error (feet) 

Landmark 

Yaw 

Lati tnde 

At Initial 

At Indicated 

Error After 

Error 



Latitude (50°) 

Final Latitude 

Correction (feet) 

(arcsec) 

(deg) 

Track 

Range 

Track 

Range 

Track 

R ange 









1 

45 


3 

0 

1. 4 

0 

1.4 

0 

10 







14 



14 



50 

1 






70 



70 



1 

40 





2. 8 



2. 8 . 



10 







28 



28 



50 

' 






140 



140 



1 

30 





5. 6 



5. 6 



10 







56 



56 



50 

< 

i 





280 



280 



1 

20 





8. 1 



8. 1 



10 







81 



81 



50 







405 



405 



1 

10 





10. 4 



10. 4 



10 







104 



104 



50 


r 





520 



520 



1 

0 





12.4 



12, 4 



10 







124 



124 



50 


1 




1 

620 

1 


620 




Notes: 1) € = 0 

2) No spacecraft librations 
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TABLE 3-4 


SENSITIVITY OF LANDMARK ERRORS TO X GYRO BIAS 
DRIFT ESTIMATE ERROR 



Notes: 1) C = 0 

2) No spacecraft librations 3 _iq 































TABLE 3-5 


SENSITIVITY OF LANDMARK ERRORS TO Y GYRO 
BIAS DRIFT ESTIMATE ERROR 


Initial 

b y 

Error 
(deg /hr) 

Final 

Latitude 

' 

(deg) 

Landmark Error (feet) 

Landmark 
Error After 
Correction (feet) 

At Initial 
Latitude (50°) 

At Indicated 
Final Latitude 

Track 

Range 

Track 

Range 

Track 

Range 

. 0015 

4 

5 

0 

0 

0 

mm 

0 

-2. 1 

. 015 









-21 



-21 

. 075 









-105 



-105 

. 0015 

40 







-4. 3 



-4. 3 

. 015 









-43 



-43 

. 075 


f 







-215 



-215 

. 0015 

30 







-8. 5 



-8. 5 

. 015 









-85 



-85 

. 075 

1 








-425 



-425 

. 0015 

2 

3 







-12.8 



-12. 8 

. 015 









-128 



-128 

. 075 









-640 



-640 

. 0015 

10 




- 



-17. 1 



-17. 1 

. 015 









-171 



-171 

. 075 









-854 



-854 

. 0015 

C 








-21. 4 



-21.4 

. 015 









-214 



-214 

. 075 

■ ' 








-1068 

i 


-1068 


Notes: 1) € - 0: 

2) No spacecraft librations 
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TABLE 3-6 


SENSITIVITY OF LANDMARK ERRORS TO Z GYRO 
BIAS DRIFT ESTIMATE ERRORS 


Initial 

b 

z 

Error 

(deg/hr) 

Final 

Latitude 

(deg) 

Landmark Error (feet) 

Landmark 
Error After 
Correction (feet) 

At Initial 
Latitude (50°) 

At Indicated 
Final Latitude 

Track 

Range 

Track 

Range 

Track 

Range 

. 0015 

i 

15 

0 

0 

0. 1 

0 

0. 1 

0 

. 015 







1. 0 



1. 0 



. 075 


r 





5 



5 



. 0015 

40 





0. 4 



0.4 



. 015 







4 



4 



. 07 5 ' 


f 





20 



20 



. 0015 

30 





1. 5 



1. 5 



. 015 







15 



15 



. 075 







75 



75 



. 0015 

2 

0 





3. 4 



3.4 



. 015 







34 



34 



. 075 







168 



168 



. 0015 

10 





5. 8 



5.8 



. 015 







58 



58 



, 07 5 

1 






292 



292 



. 0015 

0 






8. 9 



8. 9 



. 015 







89 



89 



. 075 

■ 






446 



446 




Notes: 1) e = 0 

2) No spacecraft librations 
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almost linearly with latitude change angle for angles less than 30°. 

In Table 3-5 the results for an initial error in pitch (y) gyro 
bias drift are similar to those for the roll gyro in Table 3-4 except 
that the landmark error is now in range. It is seen that the landmark 
error varies almost linearly with latitude change angle. In this case, the 
slight departure from perfect linearity is due to the inclination of the 
orbit (99°). 

In Table 3-6 the sensitivity results are shown for initial yaw (z) 
gyro bias drift error. Here it is seen that the sensitivity is much smaller 
than that for the other two gyros. 

3. 2. 1.5 Sensitivity to Attitude and Gyro Bias Drift Errors for Different 

Scan Beam Angles 

Shown in Table 3-7 is the effect of attitude and gyro bias drift 
errors on landmark error for three different scan beam angles ( e = 0, 

+ 5. 3°). The results are for a complete pass over the USA (50° to 30° 
latitude) and show the errors in the control point on the northern border 
and in an unknown landmark on the southern border. The results for 
e = o were taken from the previous tables and provide a quick comparison 
of the relative sensitivities to attitude and gyro bias drift errors. It 
should be noted that the results given for each value of e are for the case 
where both the control point and the unknown landmark were observed at 
that angle. 

In comparing the results of Table 3-7 for the three values of e 
it is seen that there is no change in the track component of the landmark 
error and very little change in the range component. Also note that the 
new range errors which appear in the right column of Table 3-7 for 
e = +5. 3° have equal but opposite polarities for the two values of e . 

Although no results were generated for the case where the scan 
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Table 3-7 


SENSITIVITY OF LANDMARK ERRORS TO ERRORS IN 
SPACECRAFT ATTITUDE AND GYRO BIAS DRIFT FOR 
DIFFERENT SCAN BEAM ANGLES FOR 20° LATITUDE CHANGE 



Landmark Error {feeP 

Landmark Error 

Initial Error 

At Start <50 u LaP 

At End (30 u Lap 

After Correction 7 feeP 


Track 

IIllSSS§!§ 

Track 

fgfgSK 

Track 

im 

For € = 0° 







1 sec in Pitch 

0 


0 

-15. 9 

0 

0 

1 sec in Roll 

-15. 9 


-14. 9 

0 


0 

1 sec in Yaw 

0 


5. 6 

0 

5. 6 

0 

. 0015 deg ^hr in 

0 


8.4 

0 

8. 4 

0 

. 0015 deg /hr in B 

y 

0 


0 

-8. 5 

0 

-8. 5 

. 0015 deg/hr in B ? 

0 


1.5 

0 

1. 5 

0 

For f- + 5. 3° 







1 sec in Pitch 

0 

-15. 9 

C 

-15. 9 

0 

0 

1 sec in Roll 

-15. 9 

0 

-14.9 


1.0 

0. 5 

1 sec in Yaw 


1.5 

5. 6 

1. 4 

5. 6 

-0. 1 

. 0015 deg 7 hr in B 


0 

8. 4 


8. 4 

-0. 1 

. 0015 deg/hr in B y 


0 


-8. 5 

0 

-8. 5 

. 0015 deg/hr in B z 


0 

1. 5 

0. 8 

1. 5 

0. 8 

For e-- - 5. 3° 







1 sec in Pitch 

0 

-15. 9 

0 

-15. 9 

0 

0 

1 sec in Roll 

(— 1 
Oi 

CD 

0 

-14.9 

-0. 5 

1. 0 

-0. 5 

1 sec in Y aw 

0 

-1. 5 

5. 6 

-1. 4 

5. 6 

0. 1 

.0015 deg /hr in B^ 

0 

0 

8. 4 

0. 1 

8. 4 

0. 1 

. 0015 deg/hr in B y 

0 

0 

0 

-8. 5 

0 

-8. 5 

. 0015 deg /hr in B z 

0 

0 

1. 5 

-0. 8 

1. 5 

-0. 8 


% aooutuc uid.i uui x euiea aaia 

V ) as control point or known landmark. 
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beam angle for the control point was different from that of the unknown 
landmark, one can see what the effect would be by analyzing the results 
in Table 3-7. For example, if e = +5. 3° for the control point and e = -5.3° 
for the unknown landmark, the initial yaw error of one arcsec would cause 
a range error after correction of -2.9 feet instead of -0. 1 feet. Note that 
the results for initial errors in pitch, roll, or gyro bias drift would be the 
same as given for e= -5.3 in Table 3-7 since these errors do not cause a 
change in the landmark error of the control point when e is changed for that 
point. 

3. 2, 1.6 Effect of Attitude Libration and Various Gyro Errors 

In Table 3-8 it seems that spacecraft attitude librations have very 
little effect on the sensitivities to attitude and gyro bias drift errors. 

Also shown in Table 3-8 are the results when gyro random drift 
and quantization noise were introduced. Although the results seem to be 
affected to some degree by gyro noise, it should be noted that the data 
happens to be that for one particular random run. A better indication of 
the effect of these noises is obtained by looking at Table 3-9 which gives 
the RMS data obtained from 20 Monte Carlo runs. Here it is seen that the 
nominal values of gyro random drift and quantization cause landmark errors 
on the order of 2-3 feet (1 <t ). 

Also shown in Table 3-8 are the effects of gyro scale factor bias 
error and input axis misalignment in the presence of spacecraft attitude 
librations. It should be noted that in the absence of attitude librations these 
errors cause only an apparent fixed bias drift which cannot be distinguished 
from true bias drift. A detailed discussion of this phenomenon can be 
found in Section 5. 5 of the Third Interim Technical Report 1 When 
attitude librations are added to the nominal attitude history, the apparent 
bias drifts vary about their mean (or fixed) values. In Table 3-9 the per- 
formance results show only the effect of the variable portion of the apparent 
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Table 3-8 


SENSITIVITY OF LANDMARK ERRORS TO ERRORS IN 
SPACECRAFT ATTITUDE AND GYRO BIAS DRIFT UNDER 
DIFFERENT CONDITIONS FOR 20° LATITUDE CHANGE 



Landmark Error (feet) 

Landmark Error 

Initial Error 


WMMMSKMl 

After Correction (feet) 


Track 


Track 

Hiisgii 

Track 

Range 








No spacecraft librations 






1 sec in Pitch 

0 

-15. 9 

0 

-15. 9 

0 

0 

1 sec in Roll 

-15. 9 

0 

-14. 9 

0 

1.0 

0 

1 sec in Yaw 

0 

0 

5. 6 

0 

5. 6 

0 

. 0015 deg/hr in B 

X 

0 

0 

8. 4 

0 

8. 4 

0 

. 0015 deg/hr in B 

0 

0 

0 

-8. 5 

0 

-8. 5 

. 0015 deg /hr in B 

z 

0 

0 

1. 5 

0 

1. 5 

0 

With spacecraft librations 






1 sec in Pitch 

0 

-15. 9 

0 

-15. 8 

0 

0. 1 

1 sec in Roll 

-15. 9 

0 

-15. 0 

0. 1 

0. 9 

0. 1 

1 sec in Yaw 

0 

0 

5. 3 


5. 3 

0. 1 

. 0015 deg/hr in B 

X 

0 

0 

8. 4 


8. 4 

0. 1 

. 0015 deg/hr in B r 

0 

0 

0 

-8. 5 

0 

-8. 5 

. 0015 deg/hr in B 

z 

0 

0 

1. 5 

0. 1 

1. 5 

0. 1 

With gyro drift noise and quantization noise, no librations 



1 sec in Pitch 

0 

-15. 9 

5. 0 

-14. 4 

3. 9 

-0. 8 

1 sec in Roll 

-15. 9 

0 

-12. 8 

0. 7 

1. 0 

1. 7 

1 sec in Yaw 

0 

0 

13. 2 

1.4 

12.3 

1. 3 

.0015 deg/hr in B 

X 

0 

0 

6. 0 

-5. 1 

5.2 

-2. 6 

. 0015 deg/hr in B 

0 

0 

3. 8 

-6. 6 

4. 0 

-9. 3 

. 0015 deg/hr in B z 

0 

0 

3. 5 

-5. 5 

3. 1 

-2.7 


Notes: (1) € = 0; (2) No gyro errors 

(31 Run with gyro noise is for one particular random run 
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Table 3-9 


SENSITIVITY OF LANDMARK ERRORS TO GYRO ERRORS AND 
NOISE FOR 20° LATITUDE CHANGE WITH ATTITUDE LIBRATIONS 


Gyro 

Error 

Landmark Error (feet) 


Landmark Error 

Error or 

Value 

At start (50° Lat) 

At end (30° Lat) 

After Correction 
(feet) 

Noise 

( la) 

■Hraraai 


Track 

m&TTFrm 

Track 


Scale 
FacLor 
Bias Error 

10 PPM 

0 

0 

0. 1 

0. 6 

0. 1 

0. 6 

Input 

Axis 

Misalign 

10 sec 

0 

0 

-4. 9 

-2. 4 

-4. 9 

-2.4 

Drift Rate 

White 

Noise 

0. 01 
deg/hr 

o ; 

0 

3.0 

2. 2 

3.0 

2.2 

Quantiza- 

tion 

Noise 

0. 10 
sec 

1. 5 

2. 0 

2.2 

2. 3 

1. 5 

1.6 

Both 

Noises 

- 

1. 5 

m 

4. 0 

2. 5 

3. 5 

2.3 


Notes: (1) Error in estimate of average bias drift rate is assumed to be zero. 

(2) Gyro errors and noise are assumed to apply equally to all 3 gyros. 

(3) Landmark error data for gyro noises represents the rms for 20 
random runs. 


(4) 


€ 


0 














































bias drift since the fixed component has the same effect as that given for 
gyro bias drift error in the preceding tables. Note also that the results 
are for only one point in the attitude libration history (i. e. t~358 seconds). 

3.2. 2 LANDMARK ERROR SENSITIVITY TO INITIAL ERRORS IN 
ORBITAL EPHEMERIS 


Tables 3-10 through 3-17 give the sensitivity of landmark errors 
to initial errors in orbital position and velocity. The magnitudes of the 
initial errors used to generate the results were 1, 3, and 6 times the 
nominal values which were 50 meters (164 feet) in down-range, 20 meters 
(66 feet) in altitude and cross-track, 0.05 meters per second (0. 164 fps) 
in altitude rate, and 0.02 meters per second (0. 066 fps) in down-range 
and cross-track rate. For brevity, the down-range and cross-track 
components of error are frequently referred to as ’range' and 'track'. 


A fairly good check can be made on the data in Tables 3-10 
through 3-16 by using the expressions for orbital position error in range 
(e RNG )# track (e TRK** and alti tude (e ALT ) in Equations 2-22 through 
2-24. Note that for a nominal attitude history and e = 0, the orbital po- 
sition errors e RNG and e TRK will also be the errors in landmark position. 
If € i 0* the altitude error will introduce a landmark error in track equal 
to e ALT sin e. 


For relatively small changes in latitude, such as across the con- 
tinental USA, the expressions in Equations 2-22 through 2-24 may be 
simplified to the following: 


e ALT " 
e RNG = 


e ALT 0 | 1+ < Wt > 2 /2 |+ 


e ALT r |COt 


+ e 


RNG 


0 


e RNG 0 \ ( ° t \ + e ALT 0 M + e RNG 0 [ 0> ^ 
1 - (*>t) |-e ALT ^ [cot j + e RNG ^ [t 


(3-6) 

(3-7) 
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e TRK 


'TRK, 


1 - (a) t) / 2 + e 


'TRK 


0 


(3-8) 


where the subscript 0 denotes the orbital ephemeris errors at the start 
of the observation pass (usually the northern USA border), t is the time 
since the start of the pass, and oo is the constant orbital rate. At a given 
point in orbit, the above errors will cause the following landmark errors: 


Down-Range Error = e RNQ (3 “9) 

Cross-Track Error = e TRK + e ALT s * n € (3-10) 

As previously indicated, the performance results in Tables 3-10 
through 3-17 are for the case where a known landmark (or control point) 
is observed at the start of the observation pass (50° latitude), and un- 
known landmarks are observed later in the pass. All of the landmarks in 
Tables 3-10 through 3-16 were observed when the scan beam angle was 
zero. In Table 3-17 data is given for the case where e 4 0. 


3. 2. 2,1 Sensitivity to Altitude Error 

In Table 3-10 the sensitivity to initial altitude error is shown for 
the case where e = 0. Note there is no error in the control point, and there 
is only an error in range for the unknown landmarks observed later in the 
pass. For small latitude changes this error varies almost linearly with 
latitude change angle and time (see first term on right of Equation 3-7). 

3. 2. 2. 2 Sensitivity to Track Error 

Shown in Table 3-11 is the sensitivity to initial error in track. 
Note that the landmark error after correction is very small in comparison 
to the error before correction for landmarks in the USA. Also note that 
the landmark error after correction varies somewhat as the square of the 
latitude change angle. 
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TABLE 3-10 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 
ALTITUDE ERROR 



Notes: 1) € = 0 

2) No spacecraft librations 
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TABLE 3-11 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 

TRACK ERROR 



Notes: 1) e - 0 

2) No spacecraft librations 
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3. 2. 2. 3 Sensitivity to Range Error 

The sensitivity of landmark errors to initial range errors is 
shown in Table 3-12. Here it is seen that the landmark error after 
correction is also small in comparison to that before correction for land- 
marks within the USA. Note that for small latitude changes, the land- 
mark error after correction varies almost as the square of the latitude 
change. 

At this point it would seem appropriate to comment on the de- 
finition of the range (or down-range) component of satellite position error 
used in this study. This error has been defined as the ’horizontal' down- 
range component of the satellite position error at a particular point in 
orbit, and not as the angular position error along the orbit (i. e. true 
anomaly error). This definition is frequently used and lends itself read- 
ily to the present treatment of satellite position and velocity in rectangular 
coordinates. Note in Equation 2-23 that for a given initial down-range 
error (e RNG ), the down-range error <e RNG ) at some other point will be 

different. However, if the initial down-range error had been defined to be 
an angular error in orbital position, the corresponding error at other points 
would have been the same and all landmarks (including the control point) 
would have been affected equally. In this case, the landmark error after 
correction would have been zero in Table 3-12. 

3. 2. 2. 4 Sensitivity to Orbital Velocity Errors 

Shown in Tables 3-13 through 3-15 are the sensitivities of land- 
mark errors to initial errors in altitude, track, and range rates. For 
small changes in latitude (i. e. with respect to 50° latitude), it is seen that 
there is less sensitivity to altitude rate error than to error in track or 
range rate. Also note that for small latitude change angles, the landmark 
error varies linearly with latitude change angle for initial errors in track 
or range rate, but varies as the square of the latitude change angle for 
initial errors in altitude rate (see also Equations 3-7 and 3-8). 
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TABLE 3-12 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 

RANGE ERROR 


Initial 

Range 

Error 

(feet) 

Final 

Latitude 

(deg) 

Landmark Error (feet) 

Landmark 
Error After 
Correction (feet) 

At Initial 
Latitude (50°) 

At Indicated 
Final Latitude 

Track 

Range 

Track 

Range 

Track 

Range 

164 

4 

5 

0 

164 

0 

163 

0 

-1 

492 





492 



488 



-4 

984 





984 



976 



-8 

164 

4 




164 



159 



-5 

492 





492 



476 



-16 

984 


r 



984 



953 



-31 

164 

30 



164 



143 



-21 

492 





492 



430 



-62 

984 

i 




984 



860 



-124 

164 

2 




164 



118 



-46 

492 





492 



354 



-138 

984 

i 




984 



708 



-276 

164 

1 




164 



84 



-80 

492 





492 



252 



-240 

984 





984 



503 



-481 

164 

C 




164 



42 



-122 

492 





492 



126 



-366 

984 

' 




984 

i 


250 



-734 


Notes: 1) * = 0 

2) No spacecraft librations 
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TABLE 3-13 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 
ALTITUDE RATE ERROR 


Initial 
Altitude 
Error 
(feet/ sec) 

.164 

.492 
. 984 

.164 

.492 

.984 

.164 

.492 

.984 

.164 
.492 
. 984 

.164 

.492 

.984 

.164 

.492 

.984 



Notes: 1) * = 0 

2) No spacecraft librations 
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TABLE 3-14 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 
TRACK RATE ERROR 



Notes: 1) * = 0 

2) No spacecraft librations 
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TABLE 3-15 


SENSITIVITY OF LANDMARK ERRORS TO INITIAL 
RANGE RATE ERROR 



Notes: 1) f = 0 

2) No spacecraft librations 
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3. 2. 2. 5 Sensitivity to Orbital Ephemeris Errors for Different Scan 

Beam Angles 

Shown in Table 3-16 is the effect of orbital ephemeris error on 
landmark error for three different scan beam angles ( £ = 0, +5. 3°). The 
results are for a complete pass over the USA (50° to 30° latitude) and show 
the landmark errors of the control point on the northern border and of the 
unknown landmark on the southern border. 

The data for each value of € in Table 3-16 is for the case where 
both the control point and the unknown landmark were observed at that 
angle. The results for e = 0 were taken from the previous tables and 
provide a quick comparison of the relative sensitivities of the different 
ephemeris errors. 

In comparing the results for the three values of e it is seen that 
a change in e affects the landmark error of the control point only in the 
presence of an altitude error. Consequently, the landmark correction, 
which uses the error of the control point, is influenced by a change in e 
only in the presence of an altitude error. For all cases, note that the 
range component of the landmark error after correction is not affected 
by a change in €. (The opposite was true in the presence of spacecraft 
attitude and gyro bias drift errors). Most of the corrected landmark 
track errors which appear when e is changed from 0° to +5. 3° are due to 
the change in £ only for the unknown landmark. 

Although no results were generated for the case where the scan 
beam angle for the control point was different from that of the unknown 
landmark, one can see what the effect would be by analyzing the results 
in Table 3-16. The only initial ephemeris error requiring consideration 
is the altitude error since the error of the control point is not affected by 
a change in e when the other ephemeris errors are present. For the case 
of initial altitude error: 
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Table 3-16 


SENSITIVITY OF LANDMARK ERRORS TO NOMINAL EPHEMERIS ERRORS 
FOR DIFFERENT SCAN BEAM ANGLES FOR 20° LATITUDE CHANGE 



Note: No spacecraft librations 
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If € = +5. 3° for the control point and € = -5,3° for the unknown landmark, 
an initial altitude error of 66 feet will cause a track error after correction 
of -12.6 feet instead of -0.4 feet. 

3*2.2. 6 Effect of Attitude Libration and Gyro Noise 

_ ' ' In Table 3-17 it is seen that spacecraft attitude librations have 
very little effect on the sensitivities to orbital ephemeris errors. Also 
shown in Table 3-17 are the results when gyro random drift and quanti- 
zation noise were present. It should be noted, however, that this is 
only for one particular random run, and that a better indication of the 
effect of gyro noise can be obtained by referring to Table 3-9 which gives 
the RMS data for 20 Monte Carlo runs. 

3. 3 EFFECT OF HIGHER HARMONICS IN EARTH'S GRAVITATIONAL 

FIELD 

In this section data is given on the effects of the higher harmonic 
coefficients of the Earth's gravitational field. The primary purpose of 
this effort was to generate data indicating the accuracy required in a model 
of the Earth's gravitational field in order to achieve a given accuracy in 
landmark location determination in the continental USA. As usual, the 
satellite orbit was initially assumed to be circular with an altitude of 540 
nautical miles and an inclination of 99 degrees. 

3.3.1 GRAVITATIONAL MODEL OF EARTH 

In Section 2. 5 some details were given on the gravitational model 
of the Earth. The reference model used in this study was the 1967 
Smithsonian gravitational model, for which the nominal values of all har- 
monic coefficients up to C g g and S g g are given in Table 3-18. The table 
also compares the 1967 values with those from other sources, including 
the more recent 1969 Smithsonian model. References 179 through 182 
were used to obtain relevant data. 
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Table 3-17 


SENSITIVITY OF LANDMARK ERRORS TO NOMINAL EPHEMERIS ERRORS 
UNDER DIFFERENT CONDITIONS FOR 20° LATITUDE CHANGE 


Initial Error 

Landmark Error (feet) 

Landmark Error 

At start (50° Lat) 

At end (30° Lat) 

After Correction (feet) 


Track 

Range 

Track 

Range 

Track 

Range 

No space craft libra 

tions 






66 ft. in Altitude 

0 

0 

0 

-24 

0 

-24 

66 ft. in Track 

66 

0 

62 

0 

-4 

0 

164 ft. in Range 

0 

164 

0 

143 

0 

-21 

. 164 fps in Alt. Rate 0 

0 

0 

-21 

0 

-21 

. 066 fps in Track Rate 0 

0 

23 

0 

23 

0 

. 066 fps in Range Rate 0 

0 

0 

22 

0 

22 

With spacecraft librations 






66 ft. in Altitude 

0 

-0. 6 

-0. 1 

-23 

-0.1 

-22. 4 

66 ft. in Track 

66 

0 

62 

0 . 1 

-4 

0. 1 

164 ft. in Range 

0 

164 

-0. 1 

144 

-0. 1 

-20 

. 164 fps in Alt. Rate 

0 

0 

-0. 1 

-20 

-0. 1 

-20 

. 066 fps in Track Rate 0 

0 

23 

0. 1 

23 

0. 1 

. 066 fps in Range Rate 0 

0 

0 

22 

0 

22 

With gyro drift noise and quantization noise, no libration 



66 ft. in Altitude 

1. 1 

2. 3 

5. 0 

-22. 9 

3. 9 

-25. 2 

66 ft. in Track 

68 

-1. 0 

64 

0. 7 

-4 

1.7 

164 ft. in Range 

0. 8 

164 

7. 6 

145 

6. 8 

-19 

. 164 fps in Alt. Rate 

0. 8 

-2. 5 

-2. 4 

-26 

-3. 2 

-23. 5 

. 066 fps in Track Rate -0. 2 

2.7 

27 

1. 9 

27. 2 

-0. 9 

. 066 fps in Range Rate 0. 4 

-2. 7 

2. 0 

16 

1. 6 

18. 7 


Notes- (II € = 0, (21 No gyro errors 

(3) Run with gyro noise is for a particular random run 
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Table 3-18 


NON- SPHERICAL HARMONIC COEFFICIENTS FOR THE GEOPOTENTIAL 

AND THEIR UNCERTAINTIES 


Coeffi- 

cient 

Coefficients of 
1967 Model 
(Normalized 
values x 10"^) 

Change 
From ’67 
to '69 
Model 

Spread in Values 
as Computed by 
Other Observers 

Assumed 

Uncertainty 

do) 

n, m 

C 

S 

(percent) 


(percent) 

(percent) 

2,0 

-1082.639 

0 

0. 001 




0. 0012 

2,2 

2. 38 

-1. 35 

1.3 




2 

3,0 

2. 565 

0 

-1. 1 




2 

3, 1 

1.71 

0.23 

16 




10 

3,2 

o 

CO 

-0. 51 

11 

-7 

to 

+ 10 

10 

3, 3 

0. 66 

1. 43 

1. 3 




3 

4,0 

1.608 

0 

-1.0 




2 

4, 1 

-0.47 

-0. 39 

18 

-8 

to 

+28 

30 

4,2 

0. 35 

0.48 

32 

0 

to 

+32 

30 

4,3 

0,92 

-0.24 

5 

-13 

to 

+ 13 

30 

4,4 

0. 04 

0. 30 

17 

-27 

to 

+77 

30 

5, 0 

0. 174 

0 

32 




30 

5, 1 

-0. 06 

i 

o 

0 

01 

41 

-38 

to 

+ 163 

40 

5,2 

0. 53 

-0.21 

23 

-25 

to 

+33 

40 

5, 3 

j 

o 

o 

0. 07 

7 

-68 

to 

+51 

40 

5,4 

-0.20 

0.02 

35 

-40 

to 

+ 185 

40 

5, 5 

0. 18 

-0. 56 

5 

-15 

to 

+ 14 

40 

6,0 

-0. 542 

0 

7 

-20 

to 

+39 

40 

6, 1 

-0. 08 

0.01 

25 

-38 

to 

+ 163 

40 

6,2 

0. 01 

-0.27 

30 

-7 

to 

+74 

40 

6, 3 

-0.04 

0.03 

. 10 

-40 

to 

+ 160 

40 

6,4 

-0. 08 

-0.48 

-17 

-23 

to 

+23 

40 

6,5 

-0.26 

-0.46 

6 

-28 

to 

+47 

40 

6,6 

-0. 02 

-0. 16 

-25 

0 

to 

+200 

40 
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The uncertainties assumed for the coefficients are given in the 
last column of Table 3-18. Due to difficulty in obtaining sufficient data to 
firmly establish the uncertainties of these coefficients, the following pro- 
cedure was used to generate the assumed values: One was to note the 
changes in the values of the C coefficients for the 1967 and 1969 Smithsonian 
models and assume that a similar change is probable in the future. The 
other was to note the spread in the values of the coefficients reported by 
other observers. The assumed uncertainties listed in Table 3-18 are pro- 
bably on the pessimistic side. Of particular interest in this table is the per- 
centage uncertainty of C 2 Q relative to that of C 2 2< This shows how well 
the most important coefficient C 9 n (frequently denoted as J_) has been 
determined. 

3.3.2 SENSITIVITY TO HARMONIC COEFFICIENTS 

Table 3-19 gives the effect on landmark location determination of 
each pair of C and S harmonic coefficients for an observational pass over the 
middle of the continental USA (50° to 30° latitude). The initial error in 
orbital position and velocity at 50° latitude was zero. The results in Table 
3-19 represent the errors which would occur in a landmark on the southern 
border if the indicated coefficient pair was not accounted for in the pro- 
pagation of satellite position and velocity. It is important to note that this 
table gives the effect on landmark location due to the coefficient pairs them- 
selves and not to errors (i.e. uncertainties) in the coefficients. It is seen 
that the landmark error due to C 2 ^ q (i. e. J 9 ) predominates by far that due 
to any other coefficient pair. The next largest error (8. 1 feet) is that due 
to the 2, 2 pair. Ihe remaining n, m pairs which cause an error of at least 
3 feet are (3,0), (3,1), (3,2), (3,3), (4,3), and (6,4). 


In Table 3-20 the landmark error is shown for various combinations 
of coefiicient pairs. The most important combination (or model) is probably 
the one containing all coefficients through C g g and S g g (except for J 2 ). It 
is seen that the omission of this coefficient combination in satellite ephemeris 
propagation will result in a landmark error of only 8. 5 feet for a complete 
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Table 3-19 


SENSITIVITY OF LANDMARK ERRORS TO INDIVIDUAL GRAVITATIONAL 
HARMONIC C AND S COEFFICIENTS FOR PASS OVER USA 


Coeffi- 

cient 

Landmark Error (feet) 

Coeffi- 

cient 

Landmark Error (feet) 


Track 

Range 

RSS 

n, m 

Track 

Range 

RSS 

2,0 

397 

-1, 503 

1, 555 

6, 2 

0 . 1 

-0. 4 

0.4 

3, 0 

-0.8 

3. 8 

3. 9 

6, 3 

0. 2 

0 

0. 2 

4, 0 

-0. 1 

1 . 1 

1. 1 

6,4 

-3. 1 

-1. 1 

3. 3 

5,0 

0 

0 

0 

6, 5 

0. 6 

-1.7 

1. 8 

6, 0 

-0. 1 

0. 5 

0. 5 

6, 6 


0.2 

0. 3 

7,0 

-0. 1 

-0.4 

0. 4 

7, 1 

-0. 3 

0. 3 

0. 4 

8,0 

0 . 0 

-0. 1 

0 . 1 

7,2 

0.7 

-1. 3 

1. 5 

9,0 

0 . 0 

0 

0 

7,3 


0. 3 

0. 5 

10, 0 

-0. 1 

0.2 

0. 2 

7,4 


0 

0. 6 

2,2 

1. 0 

8. 0 

8. 1 

7,5 

0.4 

0. 2 

0.4 

3,1 

3. 9 

-0. 3 

3. 9 

7,6 

0. 5 

-0.6 

0.8 

3,2 

1 . 0 

3. 1 

3. 3 

7,7 

0 

-0. 1 

0 . 1 

3, 3 

0. 5 

-4.4 

4. 4 

8,1 

0 

0 

0 

4,1 

-0. 8 

2. 5 

2.6 

8,2 

0 . 1 

0 

0 . 1 

4,2 

-2. 4 

-0.4 

2. 4 

8,3 

0 

-0. 1 

0 . 1 

4, 3 

-5.2 

-1. 5 

5. 4 

8,4 

-0.2 

-0.4 

0.4 

4,4 

0.7 

0. 4 

0. 8 

8, 5 

-0. 9 

-0. 1 

0. 9 

5, 1 

0 

0. 3 

0. 3 

8, 6 

-1.2 

-1.0 

1.6 

5,2 

0. 9 

-2. 1 

2.3 

8,7 

0 . 1 

-0. 1 

0 . 1 

5,3 

2.4 

0. 6 

2. 5 

8,8 

0. 2 

0 

0.2 

5,4 

-0. 5 

0. 6 

0. 8 

9, 5 

-0. 3 

0 

0.3 

5, 5 

1. 2 

-0.4 

1.3 

9, 8 

0. 3 

-0. 1 

0. 3 

6, 1 

0 . 1 

0 

0 . 1 

10,7 

-0.4 

-0. 1 

0. 4 


Note: Only the most significant coefficients among the C Q and Cm 

coefficients are listed. 9j m iUj m 
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Table 3-20 


LANDMARK ERROR VS. CHANGE IN LATITUDE FOR SEVERAL 
NON-SPHERICAL GRAVITY HARMONIC MODELS 


Lati- 

tude 

(deg) 

Lat. 

change 

Angle 

(deg) 

Landmark Error 

(feet) 

Due to all coeffs. 

thru C r 0 

6, 6 6,6 

(except Jg) 

Due to 

z, 2* °2, 2 
coefficients 
only 

Due to 
0. 001 j 2 
coefficient 
only 

50 

0 

0. 0 

0. 0 

0. 0 

40 

10 

2.4 

1. 7 

0. 4 

30 

20 

8. 5 

7. 7 

1. 5 

20 

30 

17. 8 

18. 9 

3.0 

10 

40 

31.4 

36. 2 

4. 6 

0 

50 

51. 8 

60. 3 

6. 1 

-10 

60 

81. 4 

91. 4 

7.1 

-20 

70 

121. 3 

129.7 

7.4 

-30 

80 

171. 8 

175. 1 

6. 6 

-40 

90 

233. 3 

227. 5 

4. 7 


Note: 0. 001 J 2 is approximately 100 x nominal uncertainty of J 2 
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pass over the USA. It is also seen that the coefficient pair (C 9 0 and 
S 2 2 ^ * s *h e major contributor in the previous combination. In the right 
column of Table 3-20 the landmark error is shown for 0. 001 which is 
about 100 times the nominal uncertainty of Jg. This emphasizes how 
insignificant the effect of the uncertainty is on landmark error. 

3.3.3 SENSITIVITY TO UNCERTAINTIES IN HARMONIC COEFFICIENTS 

Having presented in Section 3.3. 2 the effect of nominal C and S 
coefficients on landmark location, we now consider the effect of uncertain- 
ties in these coefficients. In Table 3-21 the effect of coefficient uncertain- 
ties is shown for various coefficient combinations. To obtain a better 
indication of performance for certain combinations, the RMS value of the 
landmark error was determined for 25 runs in which the uncertainties were 
randomized on a one-sigma basis. Most significant in Table 3-21 are the 
results obtained for all coefficients up to C g g and S g g . The RSS land- 
mark error was 1. 7 feet, which is smaller than that obtained with the Monte 
Carlo runs for gyro drift and quantization noise in Table 3-9. Also of 
interest in Table 3-21 is the fact that if the uncertainties for all C and S 
coefficients (except Jg) were 30% of the nominal values of the coefficients, 
the resulting RSS landmark error is 4. 5 feet, which is not much greater 
than 1. 7 feet. 

To summarize, it is clear from the preceding results that the 
effect of gravitational uncertainties on landmark error is insignificant. 

For this reason, the combined error studies reported on in Section 3.5 
did not include the effects of gravitational uncertainties. 

3.4 EFFECT OF ERRORS IN THE SMOOTHED ESTIMATES OF 

SPACECRAFT ATTITUDE AND GYRO BIAS DRIFT 

The data in this section attempts to show the landmark errors 
produced by the errors in some of the smoothed estimates of spacecraft 


3-35 



Table 3-21 


EFFECT OF UNCERTAINTIES IN NON-SPHERICAL GRAVITATIONAL 
HARMONIC COEFFICIENTS ON LANDMARK ERROR 


UNCERTAINTY 

in 

FOLLOWING 

Monte 

Carlo 

Run? 

Landmark Error (feet' 
Track Range RSS 

c 2 o (J 2 ) onl Y 

No 

0. 01 

0. 02 

0. 02 

C 2, 2 ; S 2, 2 ° nly 

No 

0. 02 

0. 18 

0. 18 

All coefficients 
through C 3> 3 ; S g g 

Yes 

0. 45 

0. 51 

0.7 

All coefficients 
through C 6< 6 ; S 6 _ g 

Yes 


1. 0 

1.7 

All coefficients where the 
uncertainties are assumed 
to be 30% of nominal values 
of the coefficients (J 2 not 
included) 

Yes 

2.9 

3. 5 

4. 5 


Note RMS results for 25 runs are given for all Monte Carlo runs. 
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attitude and gyro bias drift of previous SIMS studies. Presented in Table 

3-22 are many sets of errors in the smoothed estimates of spacecraft 

attitude and gyro bias drift. These are the errors which were present at 

the start of the observational pass over the USA. Most of the error sets 

1 77 

represent the errors given in the Final Report for SIMS-A configurations 

with different star mapper measurement errors and fields-of-view (FOV), 

and different data processing intervals. The error sets represent sample 

errors since they were obtained for particular sequences of random star 

mapper measurement error, random gyro drift, etc. Consequently, these 

error sets are not a statistical indication of performance of the SIMS-A 

configurations. To obtain a statistical indication of performance one must 

use the uncertainties given for the smoothed estimates in the Final 
177 

Report . There it will be found that the yaw attitude uncertainty for 
most of the SIMS-A configurations exceeds 3. 6 arcsecs (0. 001 degree), 
which was the requirement for each axis in the original SIMS studies. 

The landmark errors in Table 3-22 represent the RSS values of 
the corrected landmark errors for an unknown landmark observed on the 
southern border (30° latitude). The control point (or known landmark) 
was observed on the northern border. It is interesting to note that the 
landmark errors for the various SIMS-A configurations are almost always 
less than 10 feet. This represents one of the most significant results in 
the present effort since it indicates that for relative mapping, we are not 
as sensitive to errors in the smoothed estimates of pitch, roll, yaw, and 
gyro bias drift as was indicated for the individual errors in Section 3.2. 1. 
This is due to the existence of strong negative correlations between certain 
errors in the smoothed estimates of attitude and gyro bias drift. The 
result is a partial cancellation of the effects of certain errors on the land- 
mark location problem. Inspection of previous SIMS-A and SIMS-B data 
on the smoothed covariances of attitude and bias drift uncertainties in- 
dicate that the correlation between certain errors approaches -1 after a 
certain amount of smoothing. Evidence of such correlation can be seen in 
Table 3-22 by comparing the SIMS-A yaw errors with the corresponding 
roll bias drift errors (b x ). Note that the magnitude of the yaw error in 
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Table 3-22 


ERRORS IN SMOOTHED ESTIMATES OF ATTITUDE AND BIAS DRIFT AND 
RESULTING LANDMARK ERRORS FOR DIFFERENT SIMS CONFIGURATIONS 


Case 

Mapper 

FOV 

(deg) 

Data 

Interval 

Estimate Error in 

RSS 

Ldmk. 

Error 

After 

Correct. 

(feet) 

Attitude (arcsee) 

Bias Drift 
(10 ^ deg/hr) 

(orbits) 

Pitch 

Roll 

Y 7 aw 

b x 

b 

y 

b 

z 

SIMS- 

A ERROR 

S 








A 

4 


1 

-0.4 


4. 1 

-5. 5 

. 15 

-.60 

8. 2 

A 



2 

0. 4 


1. 6 

-1.4 

. 15 

0 

1.0 

A 



4 

-0. 3 


1.2 

-0. 9 

. 30 

15 

3. 0 

B 



1 

0. 4 

1. 8 

-15. 1 

13. 8 

-.30 

-.75 

6. 2 

C 



1 

0. 6 

-1. 5 

-7.9 

7. 1 

-.45 

. 60 

6. 1 

D 



1 

0 

-0.7 

5. 8 

-3. 9 

. 15 


10. 8 

D 



2 

0.4 

-0.5 

- 0. 8 

1. 1 

. 30 


1.7 

E 



1 

-0. 1 

0.2 

12. 3 

-12. 0 

. 45 


3. 9 

E 


► 

2 


0.7 

3. 2 

- 2. 9 

. 30 


3.2 

F 

10 

1 

1. 1 

0. 5 

7. 9 

- 8.4 

. 30 


3. 8 

F 


4 

-0. 1 

-0.7 

1. 2 

- 2. 0 

. 15 

- .30 

5. 3 

G 


4 

2 

1. 9 

-3.4 

16. 5 

-14. 7 

. 30 

. 60 

7.2 

G 


4 

6 

1.0 

1. 3 

15. 9 

-16. 7 

. 15 

-1. 35 

4. 5 

H 


1 

3. 3 

3. 1 

22. 9 

-28. 5 

.75 

-3.8 

32. 5 

H 


4 


-1. 0 

9. 2 

- 9.2 

. 15 

-0. 6 

1. 7 

SIMS- 

B TYPICAL ERRORS (see text) 






- 

- 

1 

-0. 6 

-0.4 

i. l 

-0. 5 

-. 15 

15 

3. 3 

SIMS- 

L TYPICAL ERRORS (see text) 






- 

- 

- 

2. 2 

2. 2 

23. 8 

-23.7 

.45 

29. 7 

60 


Notes: Cases B, C - different mapper and gyro noise sequences used than for A. 

Case D - no spacecraft librations; otherwise librations are present. 

Case E - large initial attitude estimation errors 

SIMS- A - Star Mapper Error (la): this is 1. 1 arcsee for Cases A 
through E. It is 2 arcsecs for Case F and 8 arcsecs for 
Cases G & H. 
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arcsecs is almost the same as the roll bias drift error in units of 10 - 
degrees per hour (or 10 arcsec per second). Also note the consistent 
reversal in polarity of these two errors. 

A fairly good check can be made on the RSS landmark errors in 
Table 3-22 as follows: Using Equations 3-1 through 3-5, the down-range and 
cross -track errors after correction can be expressed as follows for the 
case of € = 0: 


Corrected Down-Range Error = -K 60 ^ - <5^ - -KAtfib 


(3-11) 


Corrected Cross-Track Error= -K 60^ -60j 


= -K 


c(A0 ) - 1 


- s( A0 ) 


s(A0) 

fib - 

1 - c(A0) 

CO 

U U 

X 

CO 

- 


6b ^ (3-12) 
z 


The RSS of the above errors represents the desired check on the results 
in the right column of Table 3-22. 


In Equation 3-12 it is also interesting to compare the terms con- 
taining the yaw error (6ip , ) and the roll bias drift error (fib ). For the 

-3 1 x 

present orbit, co » 10 radians per second. Consequently, these two 

terms will cancel each other out if 5^ 1 (in arcsecs) = -6b (in 10" 3 arc- 

secs per second), as is almost the case for most of the SIMS-A errors 

in Table 3-22. 


The data presented in Table 3-23 shows how the landmark error 
varies for different conditions, given a particular set of errors in the 
smoothed estimates of attitude and gyro bias drift. This error set was 
taken from Table 3-22 and represents one of the worst SIMS-A configura- 
tions studied in the past. Table 3-23 shows that attitude librations do not 
have any significant effect on performance except in the presence of certain 
gyro errors. 
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Table 3-23 


EFFECT OF A SET OF SIMS-A SMOOTHED ESTIMATE ERRORS ON 
LANDMARK ERRORS UNDER DIFFERENT CONDITIONS 


Parameters 

Landmark Error (feet) 


Landmark 

Libra- 

Gyro 
Drift & 
Quant. 

Gyro Bias, 
Scale Factor 
Error & IA 

At Initial 
Latitude (50°) 

At Final 
Latitude (30°) 

Error after 
Correction (ft. ) 

tions 

? 

Noise 

? 

Misalignment 

? 

Track 

Range 

Track 

Range 

Track 

Range 

No 

No 

No 

10. 3 


21. 1 

-0. 4 

10. 7 

-0. 9 

Yes 

No 

No 

10. 3 

1; 

19. 4 

-0.4 

9. 1 

-0. 9 

No 

Yes 

No 

10. 4 

j jl&il 

23. 2 

3. 5 

13. 0 

3. 0 

No 

No 

Yes 

10. 3 


21. 1 

-0.4 

10. 7 

-0.9 

Yes 

No 

Yes 

10. 3 

I 

14. 5 

-2. 2 

4. 2 

-2.6 

No 

Yes 

Yes 

10. 4 


23. 2 

3. 5 

13. 0 

3.0 

Yes 

Yes 

Yes 

10.4 

1. 6 

16. 7 

4. 1 

6.7 

3. 9 


Notes- (1) Data above obtained with following SIMS-A error set: 

[-.03, -.7, 5.8; -.0039, .00015, .001051 arcsec, deg/hr 
These smoothed estimate errors were obtained for the midpoint 
of a 1 orbit data interval with 1. 1 arcsec mapper sensor error and 
4° optics FOV. 

(2) Scan beam oriented downwards, € = 0. 

(3) No initial navigation errors. 

(4) When gyro noise is present, data given is for the rms of 
25 Monte Carlo runs. 
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3. 5 COMBINED EFFECT OF ERRORS IN SPACECRAFT ATTITUDE 

AND ORBITAL EPHEMERIS 

Until now, consideration has only been given to the separate 
effects of errors in spacecraft attitude and orbital ephemeris. In this 
section, these errors are combined to determine the overall effect on 
landmark error. No consideration was given to gravitational coefficient 
uncertainties since their effects were found to be relatively small for 
observational passes over the continental USA. Although the data pre- 
sented in this section may not give as much insight into the effects of in- 
dividual error sources, it has been generated as a matter of interest. 

To implement this phase of the study, the following errors were 
adopted as nominal for spacecraft attitude and gyro bias drift at the start 
of the observational pass: 

[2, 2, 4; -.0030, -.0015, .0015 ] arcsec, deg/hr 

which are the respective errors in pitch, roll, yaw, and the bias drift 
of the X, Y, and Z gyros. The 4 arcsec error in yaw approximately 
corresponds to the 0. 001 degree error originally specified for the SIMS 
trade study. The reason for the value of 2 arcsecs in pitch and roll is 
simply a recognition of the better performance observed for these two 
angles in the previous SIMS studies. The gyro bias drift errors were 
chosen to be compatible with the attitude errors in accordance with the 
smoothed results of previous studies. 

The nominal errors adopted for the orbital ephemeris were the 
same as those given in Section 2.4. As previously indicated in Section 2.4, 
the data used to establish these values was somewhat limited but was con- 
sidered to be sufficient for most of the studies. However, in the present 
case, it would have been nice to have had one or more sets of errors which 
are representative of the present ground tracking and batch -pro cessing 
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techniques. It is felt that this data would reflect the correlation between 
certain ephemeris errors, just as was the case for the smoothed estimates 
of spacecraft attitude and gyro bias drift in the SIMS studies. It would also 
have been desirable to have one or more covariance matrices of the 
ephemeris errors so that randomized sets of errors could be generated. 

In the present effort, the nominal values selected for the ephemeris errors 
were: 


[ -20, 20, 50; -.05, . 02, .02 ] meters, meters/sec 

which are respectively the errors in altitude, track, range, altitude rate, 
track rate, and range rate. The relative magnitudes and polarities of 
certain errors were adopted to reflect certain error relationships which 
have been encountered in other studies at CSDL. 

3.5.1 COMBINED EFFECT OF NOMINAL ATTITUDE AND EPHEMERIS 

ERRORS 

Tables 3-24 through 3-27 and Figures 3-1 through 3-5 present data 
on the separate and combined effects of nominal errors in the ephemeris 
and of nominal errors in spacecraft attitude and gyro bias drift. Table 
3-24 gives data on the uncorrected landmark errors versus latitude for 
an observational pass over the USA, This data is illustrated graphically 
in Figures 3-1 and 3-2. It is seen that the nominal ephemeris errors 
make a greater contribution to the uncorrected landmark errors than the 
nominal errors of attitude and gyro bias drift errors. Shown also is the 
effect of different scan beam angles. Note in Table 3-24 and Figure 3-1 
that a variation of the scan beam angle, c , has almost no effect on the 
down-range component of landmark error when only ephemeris errors 
are present, and that the same is true for the cross-track component of 
landmark error when only attitude and gyro bias drift errors are present. 
Thus when we inspect the results for the combined error sources in 
Figure 3-2, we know from the above that the spread in the down-range 
error curves for different scan beam angles is due to the attitude and gyro 
bias drift errors; while the spread in the cross -track error curves is due 
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Table 3-24 


UNCORRECTED LANDMARK ERRORS VS. LATITUDE DUE TO THE PRINCIPAL 
NOMINAL ERROR SOURCES FOR DIFFERENT SCAN BEAM ANGLES 


Latitude 

Landmark Errors in Feet Due to Nominal 
Initial Errors In 



Ephemeris 

Attitude and 
Bias Drift 

Both Cases 
Combined 

(deg) 

Track 


Track 


Track 


For € = 0 







50 

65. 6 

164. 6 

-31. 8 

-31. 8 

34. 9 

135. 1 

45 

71. 7 

175. 9 

-30. 6 

-34. 4 

43. 9 

142.4 

40 

76. 9 

187. 1 

-29. 1 

-36. 7 

48. 7 

151. 5 

35 

81. 2 

198.2 

-27. 1 

-38. 7 

55. 1 

159.9 

30 

84. 6 

209. 4 

-24. 2 

-40. 4 

60. 7 

168.4 

For e = +5. 3° 







50 

59. 6 

164. 6 

- -32.0 

-25. 9 

28. 7 

141. 1 

45 

65. 7 

176.0 

-30. 8 

-28. 0 

37. 7 

149. 1 

40 

71. 0 

187. 1 

-29. 3 

-30. 0 

42. 6 

158. 9 

35 

75. 4 

198. 2 

-27. 3 

-31. 5 

49. 1 


30 

79. 0 


-24. 5 

-32. 8 

54. 8 

176. 9 

For e = -5. 3° 







50 

71. 8 

164. 5 

-32.2 

-37. 7 

40. 7 

129. 1 

45 

77. 8 

175. 9 

-31. 0 

-40. 7 

49. 7 

135. 6 

40 

83. 0 

187. 0 

-29. 5 

-43. 5 

54. 5 

144. 2 

35 

87. 2 

198. 2 

-27. 3 

-45. 9 

60. 9 

151. 9 

30 

90. 4 

209.4 

-24. 3 

-48. 0 

66. 4 

159. 9 


Note* (D Runs made with spacecraft librations 

(2) Data for combined errors includes effect of gyro errors and noise 
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, — CROSS TRACK ERROR J 

BIAS DRIFT ERRORS 


LATITUDE (DEGREES) 


Figure 3-1 Uncorrected Landmark Errors Versus Latitude 
Due to Nominal Initial Errors in Ephemeris, 
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2 Uncorrected Landmark Errors Versus Latitude 
Due to Nominal Initial Errors in Ephemeris, 
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to the ephemeris errors. This result could be of significant value in 
future studies involving the use of known landmarks to estimate both the 
spacecraft attitude and orbital ephemeris. 

In Table 3-25 the corrected landmark errors are presented for 
the combined ephemeris and attitude error cases of Table 3-24. For the 
first time in this report, data is presented showing the effect of using a 
control point (or known landmark) observed at a latitude other than 50°. 

In Table 3-25 the corrected landmark error is given for an unknown land- 
mark using the correction data obtained for a control point either at 50°, 
40°, or 30° latitude. It is clear from Table 3-25 that the use of a control 
point at 40° latitude (midpoint of pass) results in a maximum landmark 
error which is about one half of that obtained when the control point is at 

50° or 30° latitude. Graphical illustration of the results for a control point 
at 40° latitude is shown in Figure 3-3. It should be noted in Figure 
3-3 that the polarity of the cross-track error curves was reversed 
since they would have coincided closely with the down-range error curves. 


3.5.2 COMBINED EFFECT OF NOMINAL EPHEMERIS ERRORS 

AND 3X NOMINAL ATTITUDE ERRORS 

As a matter of interest, data was also generated showing the 
effect of combining nominal ephemeris errors with attitude and gyro bias 
drift errors which were three times nominal. The data for this case is 
shown in Table 3-26 and Figures 3-4 and 3-5. Of particular interest is 
the increase in the spread of the down-range (or range) error curves for 
different scan beam angles (see Figures 3-2 and 3-5). From the preceding 
discussion we know that the spread in the down-range error curves is due 
to the attitude and bias drift errors. In particular, it is due to the presence 
of a yaw attitude error at the time of observation of the unknown landmark. 
The yaw error did not affect the control point since e = 0. 


3-46 



Table 3-25 


LANDMARK ERRORS VS. LATITUDE FOR NOMINAL INITIAL 
ERRORS IN EPHEMERIS, ATTITUDE AND BIAS DRIFT 


Indicated 

Latitude 

(deg) 

Landmark Error (feet) 

At Indie 
Latitu 

ated 

de 

After Correction Using Known Landmark 
Along e = 0° at 


50° Lat. 

40° Lat. 

30 

msSml 

Track 


Track 

ESSHI 

Track 


Track 

BigffBBl 





mmmammmm 





ii 

o 

o 









50 

34.9 

135. 1 

0 

0 

-13. 8 

-16. 5 

-25. 8 

-33. 3 

45 

43. 9 

142.4 

9.0 

7. 3 

-4. 8 

-9. 2 

-16. 8 

-26. 0 

40 

48.7 

151. 5 

13. 8 

16. 5 

0 

0 

-12. 0 

-16. 9 

35 

55. 1 

159.9 

20.2 

24. 9 

6.4 

8.4 

-5.6 

-8. 5 

30 

60.7 

168.4 

25.8 

33. 3 

12. 0 

16. 9 

0 

0 

e=+ 5. 3° 









50 

28.7 

141. 1 

-6. 2 

6. 0 

-20.0 

-10. 5 

-32. 0 

-27. 3 

45 

37.7 

149. 1 

2. 9 

14. 1 

-10. 9 

-2.4 

-23. 0 

-19. 3 

40 

42.6 

158.9 


23. 8 

-6. 1 

7. 3 

-18. 1 

-9. 5 

35 

49. 1 

168.0 


32. 9 

0.4 

16. 5 

-11. 6 

-0.4 

30 

54. 8 

176. 9 


41. 8 

6. 1 

25.4 

-5.9 

8. 5 

C 9° 

c = - 5. 3 









50 

40.7 

129. 1 

5. 8 

-6. 0 

-8. 0 

-22. 5 

-20.0 

-39. 3 

45 

49.7 

135. 6 

14. 8 

0. 5 

1 . 0 

-15. 9 

-11. 0 

-32. 8 

40 

54. 5 

144. 2 

19.6 

9. 1 

5.8 

-7. 3 

-6. 2 

-24. 2 

35 

60. 9 

151. 9 

26. 0 

16. 8 

12.2 

0. 3 

0. 2 

-16. 5 

30 

66. 4 

159. 9 

31.5 

24. 8 

17. 7 

8.4 

5. 7 

-8. 5 
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LANDMARK ERROR (FEET) LANDMARK ERROR (FEET) 



LATITUDE (DEGREES) 

Figure 3-3 Corrected Landmark Errors Versus Latitude 

for Nominal Initial Errors in Ephemeris, Attitude 
and Bias Drift. 



Figure 3-4 Corrected Landmark Errors Versus Latitude 
for Nominal Initial Ephemeris Errors and 
3X Nominal Attitude and Bias Drift Errors. 
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Table 3-2 6 


LANDMARK ERRORS VS. LATITUDE FOR NOMINAL INITIAL 
EPHEMERIS ERRORS AND 3X NOMINAL ATTITUDE AND BIAS 

DRIFT ERRORS 


Indicated 

Latitude 

(deg) 

Landmark Error - feet 

At Indicated 
Latitude 

After Correction Using Known 
Landmark along c = 0° at 

50° Lat. 

40° Lat. 

30° Lat. 

Track 


Track 


Track 


Track 

n 

€ = 0 









50 

CO 

CO 

CM 

1 

71.4 

0 

0 

-19. 2 

- 6.6 

-41. 0 

-16.2 

45 

-IT. 3 

73. 6 

11. 4 

2. 2 

- 7.8 

- 4. 5 

-29. 6 

-14. 1 

40 

- 9.6 

78. 1 

19.2 

6. 6 

0 

0 

-21. 8 

- 9.6 

35 

1.0 

82. 5 

29.7 

11. 1 

10. 5 

4.4 

-11. 3 

- 5.2 

30 

12. 3 

87.7 

41. 0 

16. 2 

21.8 

9.6 

0 

0 

{:= +5.3 0 









50 

-35.4 

89.2 

-6.6 

17. 8 

-25. 8 

11.2 

-47. 6 

1.6 

45 

-23. 8 

93. 1 

5.0 

21.6 

-14.2 

15. 0 

-36. 0 

5.4 

40 

-16. 0 

98. 9 

12.8 

27. 5 

- 6.4 

20. 9 

-28. 2 

11.2 

35 

- 5. 5 

104. 9 

23.2 

33. 5 

4. 1 

26. 9 

-17.8 

17.2 

30 

5. 8 

111. 4 

34. 5 

39. 9 

15. 3 

33. 3 

-6.5 

23.7 

e = -5.3° 









50 

-23.6 

53. 6 

5. 1 

-17.8 

-14. 1 

-24.4 

-35. 9 

-34.0 

45 

-12. 3 

54. 1 

16. 5 

-17. 3 

- 2.7 

-23.9 

-24. 6 

-33. 5 

40 

- 4.4 

57.2 

24. 3 

-14. 3 

5. 1 

-20. 9 

-16.7 

-30. 5 

35 

6. 3 

60. 0 

35. 0 

-11.4 

15. 8 

-18. 0 

- 6.0 

-27. 6 

30 

17.7 

64. 0 

46. 5 

- 7.5 

27. 3 

-14. 1 

5. 5 

-23.7 
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Figure 3-5 Uncorrected Landmark Errors Versus Latitude 
for Nominal Initial Ephemeris Errors and 3X 
Nominal Attitude and Bias Drift Errors. 
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It is important to note that Figures 3-3 and 3-4 can be used to 
obtain similar data for the case where the control point is located at a 
latitude other than 40° but still observed at e= 0. For example, to 
obtain the results for the case where the control point is at 50° latitude, 
one simply translates vertically the curves in the figure until the down- 
range and cross-track error curves for c= 0 pass through zero for 50° 
latitude. 


Table 3-27 gives data on uncorrected and corrected landmark 
errors for various error configurations. The important things to note 
here are: (1) The landmark errors after correction are three to five 
times greater with nominal ephemeris errors than with nominal attitude 
errors, when these are considered separately. (2) The relative effect 
on landmark errors after correction is much less when the nominal 
attitude errors are increased by a factor of three then when the nominal 
ephemeris errors are increased by the same factor. 
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Table 3-27 


UNCORRECTED AND CORRECTED LANDMARK ERRORS FOR VARIOUS 

ERROR CONFIGURATIONS 


Nominal Errors for 

Landmark Error (feet) 

Landmark 
Error After 
Correction (ft! 

Ephem- 

eris 

Attitude 
& Bias 
Drift 

Gyro 
SFE & 
IA Mlm. 

Gyro 
Drift & 
Quant. 
Noise 

At Start 
(Lat= 50°) 

At End 
(Lat= 30°) 

Track 

Range 

Track 

Range 

Track 

Range 











IX 

0 

0 

0 

65. 6 

164. 6 

84. 6 

209. 4 

19. 0 

44. 8 

IX 

0 

IX 

IX 

66.7 

166. 9 

84. 9 

208. 8 

18.2 

41. 9 

0 

IX 

0 

0 

-31. 8 

-31.8 

-24. 2 

-40. 4 

7. 6 

-8. 5 

0 

IX 

IX 

IX 

-30.7 

-29. 5 

-23. 8 

-41. 0 

6. 9 

-11. 4 

0 

0 

IX 

IX 

1 . 1 

2. 3 

0. 4 

-0. 6 

-0. 7 

-2. 9 

IX 

IX 

IX 

IX 

34. 9 

135. 1 

60. 7 

168.4 

25. 8 

33. 3 

IX 

3X 

IX 

IX 

-28. 8 

71. 4 

12. 3 

87.7 

41.0 

16.2 

3X 

IX 

IX 

IX 

166. 1 

464. 3 

229. 8 

587. 1 

63.7 

122.9 

3X 

3X 

IX 

IX 

102.4 

400. 6 

181. 4 

506. 4 

78. 9 

105. 8 

Following runs are 

without librations 







IX 

IX 

IX 

IX 

34. 9 

134. 5 

67. 0 

170. 7 

32. 1 

36. 3 

IX 

0 

0 

0 

65. 6 

164. 0 

84. 4 

210. 0 

18. 8 

46. 0 

0 

IX 

0 

0 

-31. 8 

-31. 8 

-22. 8 

-40. 4 

9. 0 

-8. 5 


Note (11: 
( 2 ): 

'3) 

(4) 

(5) 


3X denotes that nominal errors are multiplied by a 3 times factor, etc. 

Correction is made assuming known landmark (or control point! 
is located at 50° latitude and along e =0. 

Runs with spacecraft librations unless otherwise noted. 

e = 0. 

Runs with gyro noise made for a particular random run. 
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SECTION 4 


CONCLUSIONS 

The primary objective of the landmark location study was to 
determine the effects of errors in spacecraft orbital ephemeris, space- 
craft attitude, and gyro bias drift on the ability to determine the locations 
of unknown landmarks with respect to known landmarks (or control points) 
in the payload sensor imagery. Most of the data was generated for obser- 
vational passes over the continental USA. It is felt that sufficient data 
was generated to clearly show the sensitivity of landmark location deter- 
mination to errors in each component of the orbital ephemeris, spacecraft 
attitude, and gyro bias drift. 

In the present study it was assumed that only one known landmark 
(or control point) was observed during a pass over the USA. The data in 
this report shows what the effects of ephemeris and attitude errors are on 
determining the locations of other arbitrary features (observed during the 
pass) with respect to this control point. In most cases the control point 
was arbitrarily selected to be on the northern USA border in order to show 
the maximum error in determining the relative location of an unknown 
feature (or landmark) on the southern border. It is important to note that 
errors of about the same magnitude would occur if the control point was on 
the southern border and the unknown landmark was on the northern border. 
As shown in Section 3. 5. 2, the error for an unknown landmark on the 
northern or southern border would have been about half as large if the 
control point had been located halfway across the USA (i, e. at 40° latitude). 
It is obvious that the relative landmark errors can be reduced by using 
more than one control point in an observational pass. The data in Section 3 
provides some indication of how close the control point must be to the un- 
known feature in order to achieve a certain minimum accuracy in relative 
mapping. However, it should be pointed out that some of the data in Section 
3 does not give a statistical indication of performance since the results are 
for discrete sets of errors in orbital ephemeris, satellite attitude, and gyro 
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bias drift. A complete statistical analysis would probably be appropriate 
in the future after better statistical data has been obtained for the orbital 
ephemeris errors (including correlation between the errors), and possibly 
after a particular SIMS configuration has been identified for this type of 
mission. 


One surprising result of this study is the apparent decrease in 
the effect of certain attitude and bias drift errors on relative landmark 
location determination when these errors happen to be the errors in the 
smoothed estimates of attitude and bias drift. By using smoothed (or 
filtered) estimates of attitude and gyro bias drift, one takes advantage of 
certain relationships established between the errors in these estimates 
which tend to decrease the effect of certain errors on relative landmark 
location determination. This was clearly shown to be the case for errors 
in yaw angle and roll bias drift for SIMS-A in Section 3.4. There it was 
shown that these two errors tend to cancel each other out. At this point it 
would seem appropriate to review the history associated with the yaw 
attitude error in the overall SIMS trade study. At the beginning of this 
study the desired accuracy in estimating pitch, roll, and yaw was 0. 001° 
(la) each, although it was felt by many that the yaw accuracy did not have 
to be as good as that for pitch and roll since yaw error had very little effect 
on landmarks observed near local vertical. In the SIMS-A studies 169 * 177 
it was found that the desired yaw accuracy was also the most difficult to 
achieve because of the nature of that system. In Section 3. 2 it is shown 
that the yaw error indeed does not introduce as much landmark error as 
pitch and roll when the position of that landmark is being determined solely 
with satellite position and attitude data. However, when one attempts to 
determine the location of an unknown landmark with respect to a known 
landmark, the landmark error is more sensitive to yaw error than it is to 
pitch or roll error. Note that these are the results obtained for the in- 
dividual errors in pitch, roll, and yaw. In looking at the individual sen- 
sitivity data of Section 3.2, it would therefore seem that the relatively 
larger yaw errors of SIMS-A would make a significant contribution to the 
relative landmark error. Yet, this is shown in Section 3.4 to not be the 
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case when the attitude and bias drift errors are those obtained as the result 
of smoothing. This conclusion would seem to be a very good argument 
tor having a SIMS onboard the satellite for this type of mission, since this 
would permit smoothed estimates of attitude and bias drift to be obtained. 
Although no smoothed data was available for ephemeris errors, it would 
also be interesting to see what effect such data would have on landmark 
errors. It is possible that smoothed estimates of satellite position and 
velocity may have an effect similar to that obtained for smoothed estimates 
of attitude and gyro bias drift. 

Another important result of this study is that obtained on the 
effect of higher harmonics in the Earth’s gravitational field. For land- 
mark location determination in a single pass over the continental USA, it 
is shown in Section 3. 3 that a maximum landmark error of about 10 feet 
will occur when using a gravitational model consisting of the central force 
field and only the well known harmonic J 2 . The above error is that 
for an unknown landmark on the southern or northern border when the 
known landmark (or control point) is on the opposite border. If the control 
point happens to be mid-way across the USA (40° latitude), the above error 
will be half as large. Sufficient data is presented in Section 3. 3 to show 
the separate effects of additional harmonics when greater accuracy is de- 
sired. The results in Section 3. 3 show that the effect of the uncertainties 
in the harmonic coefficients is much smaller than that of the harmonic 
coefficients themselves. 
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SECTION 5 


RECOMMENDED FUTURE STUDIES 


The present landmark location study was primarily concerned 
with determining the effects of spacecraft orbital ephemeris and attitude 
errors on the ability to determine the locations of unknown landmarks with 
respect to known landmarks when a relatively simple but effective tech- 
nique of landmark correction is used.' Although the results for each case 
in this report are for just one known landmark per observational pass, it 
is felt that sufficient data was generated to show the effects of known land- 
mark location at different points in the imagery. No attempt was made to 
use an optimal technique of correction for relative landmark location since 
this was beyond the scope of this study and was considered to be more 
appropriate in future studies. 

Although sufficient sensitivity data was generated in the present 
study to enable one to estimate what the performance would be for various 
combinations of errors in orbital ephemeris and satellite attitude, it is 
possible that additional data may be desired in the future using errors 
which are more realistic than those defined to be nominal in the present 
effort. This would, of course, depend on decisions as to what SIMS con- 
figuration is adopted. Also, as indicated in Section 4, it would be in- 
teresting to see what the effect would be on landmark location determination 
when smoothed estimates of orbital ephemeris are used. 

One area which is strongly suggested in future studies is the use 
of known landmarks in the imagery to update the orbital position and velocity, 
the satellite attitude, and the gyro bias drift. This represents an update 
of twelve parameters and should be done with the Kalman filter and/or an 
optimal smoother. It is felt that this approach may indirectly provide the 
optimal technique of landmark correction mentioned earlier. 
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Such a study should initially be a covariance analysis with consideration 
being given to various gyro package configurations, landmark distributions, 
etc. Consideration should also be given to the addition of star measurements 
by various candidate star sensors. Eventually, a full state simulation 
should be performed on the most promising configuration (or configura- 
tions) with optimal estimates being made of both the state and the co- 
variance matrix of the state estimation error. 

In addition to the above, the Active and Passive Landmark 
Mechanization Study, which was not performed in the present study effort, 
should be done in the near future because of its possible value to automatic 
ground data processing of thematic mapper imagery. 
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APPENDIX A 


Report 

Page 

2-4 


2-16 

2-17 


ERRATA FOR FINAL REPORT R-741 


The following corrections should be made in the Final 
77 

of which this report is an addendum: 

Correction 

The value of given at top of page should 
be 0. 478 instead of 0. 005. Also remove 
1 -3cr 

In Equation 2-3 draw line under aVp to denote 
vector. 

Place dot above P on left side of Equation 
2-32. 


Place dot above U on left side of Equation 
2-36. 
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